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ABSTRACT 


Large Reynolds number approx im.at ions are examined 
for certain internal separated flows hy considering various 
limits. It is argued that the equation governing the flows 
except possibly for some subregions, has the form of the 
boundary-layer equation, as the inertial and viscous forces 
are of the same order. A method of calculation is then 
developed which uses the eigenfunctions of the linear 
problem governing the far-d own stream flow. It leads to a 
set of coupled, first order, ordinary differential equations 
It is first tested for two cases for which results ba.sed cn 
the Navier-Stokes equations are available. Details of the 
flow including streamlines in the recirculating regions are 
well predicted. It is concluded that the equation represent 
the large Reynolds number behaviour of the flovis considered. 
Plow problems with three types of wall geometry are then 
solved. They are a channel having symmetric sudden expansio 
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a channel having asjTiimetric sudden expansion and a "base 
placed symmetrically in a channel. Calculations are 
performed for various values of the ratio of the upstream 
channel width to that downstream. Also two types of entr3/' 
conditions representing a long and a short inlet are consi- 
dered. Eesults which are presented include centreline 
velocity, velocity profiles, development length, streamlines, 
skin friction, pressure distribution and pressure recovery 


coefficient . 



CHAEEEE 1 


lETRODUCTIOlT 

Separated flows constitute an important tut difficult 
area of fluid mechanics. Bound ary- layer approximations 
provide substantial simplifications in attached flows at 
large Reynolds number. In particular, the outer flow can be 
obtained to the lowest order, by an inviscid calculation in 
which the condition of zero normal velocity at the wall is 
satisfied. This is because the bound ary- layer and the wahe 
where the viscous forces and the inertial forces are of the 
same order, are' rather thin and their effect on the outer 
flow is of a higher order. Separated flows, on the other, 
hand, have quite thick rotational regions and the outer flow 
is affected substantially. Consequently, the outer flow 
calculation presents certain difficulties which have yet to 
be fully resolved. 

The present thesis deals Vi/ith a simpler clans of 
separated flows v^^hich has one ingredient in common w'ith 
external separated flows, namely that they have thick 
rotational regions. We consider the flows confined between 
t?/o walls -(Bigure 1), v'^hich belong to the category of 
internal separated flow?s. One such flow is the flovj in a 
channel having a sudden expansion in which recirculating 
regions are formed downstream of the expansion. As we shall; 
see in Chapter 2, there are two types of internal separated : 
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flows which differ in their streamwise length scales. 


limiting equations for various longitudinal and 
transverse length scales are, examined as R (Cnapter 2). 

It turns out that for the type of in'cernal separated flow 
problem referred to earlier, the streamwise length scale is 
of order R. As the inertial forces are of the order of the 
viscous forces in the limit considered, the equation has the 
form of the -hound ary-layer equation. However, the longi- 
tudinal and transverse length scales of the present problem 
ajce different from the classical boundary-layer scales. 

Also, the pressuT'e gradient is an additional ujiloiown which 
is determined by using four boundary conditions arising from 
the v>;alls unlike the boundary-layer problem. 

A oalculation procedure of the equation, which is 
developed in Chapter 3, ex;ploits the ultimate development of 
the flow into Poiseuille flow by taking the solution in 

the form 




(x,y;R) = '-kCy) + a (x/E) dl; (y) 

1 • 

1 


YJhere x and y are the strearm^/ise and transverse coordinates 
and is the stream function. Usually in an integral 
relations method, the expansion set din and the weighting 
functions are chosen for their anal37tical properties and 
computational simplicity. Y'e however, adopt .a. different 
approach here. The flow at large downstream distances is 
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a small of the Poiseuille flow vJhich leads to 

an eigenvalue p^^oblem. Its eigenfunctions form an efficient 
expansion set as the resulting set of coupled ordinary 
differential einations tend to get decoupled with increasing 
X. The -weighting functions are the natural ones arising 
from the eigenvalue problem. The given flow problem then 
reduces to integrating these differential equations which 
are coupled quasi-linear and of first order. The initial 
condition has obtained by matching the far -flow field 

with the ne ar—f 1°''''' field. In the absence of any rigorous 
description of f^ie near-flow field, a hypothesis is made to 
give the initial condition. 

The schenie of calculation is first applied to two test 
cases for which numerical results based on the complete 
Havier-Stoloes equations are available, and can be compared 
with the preseJ^f calculation. The first problem to be taker 
up is the now classical entry flow problem (Chapter 4-). Ite 
simple geometri’' forms a suitable testing ground for theories 
and methods amd hence is considered here, although no 
recirculating regions are present. The second test case is 
the flow in a, channel having a symmetric sudden expansion 
w;ith an expan^f®^^ ratio of two. Comparison of the present 
calculation wfth that based on the havier-Stoke s equations 
for large E (Snng and Macagno , 1966) is found to be very | 

good . Not onXy ftie gross quantities of interest such as thi 
length of tiie recirculating region and the total recirculat; 
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are foimd to agree well, but also the streamline pattern is 
well predicted. These results are discussed in Chapter 5. 

The above comparison shovvs that the present method of 
expansion into eigenfimctions of the far downstream flow 
provides a suitable calculation procedure , and it also 
demonstrates the applicability of the limiting equation. 

We thus conclude that the streamwise length scale for such 
flows is of order R, the viscous term and the inertia term 
are of the same order, and the governing equation is 
parabolic despite the presence of recirculating regions. 

Flow problems viith three different types of wall 
geometry are solved by using the present method (Cee 
Figure 1). One type illustrates the applicability of the 
method to asymmetric flow problems. Another type models a 
base between two walls. Calculations are made for various 
values of the geometrical parameter h, and are included in 
Chapter 5. The results of these calculations are believed 
to be new. 

As the method uses an eigenfunction expansion, the 
calculation, in principle, should be continued with increas- 
ing number of terms until a satisfactory convergent trend is 
obtained. However, we have used a rather small inumber of 
termis in the calculations. For example, the favourable 
comx'arison in the cane of the sjnmmetric sudden 



expansion mentioned eanlier is obtained v>/ith only three 
terms. The maximum number of terois used in the present 
calculations is therefore limited to five. Obviously the 
choice of the number of terms v^ould depend on the v\/all 
geometry, the entry profile and the desired accuracy. 

If the number of terms to be used in a calculation is 
fixed, a difficulty may be encountered in the integration 
of the ordinary differential equations. 'This difficulty is 
essentially due to the non-linear evolution term and is 
a general feature vlienever a moment method or an integral 
rela-tions method is applied to the present type of partial 
differential equation. This issue is discussed in Section 
3.5. Phase plane analysis of the ord inary d ifferent ial 
equation when two terms are used , is given to illustrate 
this feature. 

We have restricted our attention to steady laminar 
two-dimensional flow of a Hewtonian incompressible fluid, 
although the basic elements of the analysis would be appli- 
cable to a larger class of problems. We have however chosen 
to apply the method to three tjj^pes of v^all geometry, for a 
range of the parameter h, two tr^'pe of entry conditions, and 
three choices of the number of terms with a view to obtain 
insight into the performance and the applicability of the 
method . 
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LIMIT EQUITIOIT FOR LARGE REIROLDS NUIBER 

2.1 GOYERHIRG EQUATION 

The equation governing the two-dimensional plane 
steady laminar motion of a Newtonian incompress idle fluid 
is 

J( v}. ) = R'V^V'^ , (2.1) 

Yihere is the stream function. The two cartesian coordi- 

2 

nates are denoted by x ' and y, and y and J denote the 
laplacian and the Jacobian operators respectively. The 
above eque.tion is written in the non-dimensional form, R 
being the Reynolds number. Equation (2.1), known as the 
vorticity transport equation, is obtained from the two 
momentum equations and the equation of continuity. The 
ahove equation, without any approximation , will be denotcsd 
by USE . 

2.2 LilRGE REYNOIDS NUlvIBER LIMITS 

Let the scales of x, y and y be denoted by L , L 
and if respectively. Then the orders of the various terms 
in equation (2.1) are 

( ( Y/rl^, t'/RL^Ly, . 

The first bracket shows the orders of inertia terms, while 
the orders of viscous terms are shown in the second bracket. 
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Terms of the lowest order in either group will depend on 

the order relation "between the two length scales 1 and L . 

Therefore three sets of limits are possible depending on 

whether r (' j L , L 'C: I . Pcssihle limit 

X j y X X j 

equations are given "below and each limit equation has "been 
designated "by E, Y, B or IT stand in.g for inviscidy viscous, 
"bound ary- layer or Navier-Stokes types of limits. 


E 


IT 


V 


E 

y 


B 

y 


J(c/ ^ 'p , V ) = 0 

J(\7^ , y ) = R""' 77^ \7^ ’4^' 

i > 

? p 

X? ^ ^ = 0 


L C. I 

y X 


J( '4'' ' , = 0 


J( vi -' , 1' ) 

' yy * ^ 


P 'sj> 


’ jjjj 



' jjyy 


0 


The order sjmihols used here have the following meaning. 
1 'X-' J-i 5 0 -illim 1 /li ^ OQ 

^ 1 H-co ^ y 

I . Cl 1^^ : lim I /I = 0 

- R-^co 7 


X 


y 
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Ii c: 1 
X V 


J( ) = 0 

M XX i ^ 

j( '..U , ) = R-^ 

^ I XX ’ ^ I xx> 3 ^ 


xxxx 


( 2 . 2 ) 


Suffixes X and j following ^ denote respective partial 
derivat Ives . 


limit equations E and Y can le obtained from B,., 

y y y 

and E and 7, from B . limit equations E and E are 

X X X y 

contained in E and V contains and Y^. Distinguished 
limits are those corresponding to E, Y, B and B and E is 


the complete equation of motion. 


E is the Euler equation governing inviscid fluid flow. 
Stokes equation in which only viscous terms are retained is 
designated by Y. Interestingly, Stokes equation is obtained 
for large Reynolds number R. However, the appropriate 
Reynolds number R goes to zero. B,^ is the familiar 

t) 

bound ary- layer equation differentiated once with respect to 
j to eliminate the pressure term. E is sometimes referred 
to as the inviscid bound ary- layer equation. 3 and E, are 

X X 

X counterparts of B and E respectively. 


The domains of validity of the various limit equations 
can be readily visualised v;ith the help of a three-dimen- 
sional plot with axes 1, , 1 and ")f (Figure 2). This plot 

X y 
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is of a special natirce as the axes are order axes. 'That is, 
a roirit, say (Pl,a/E, 1), denotes LaJ R, I . e VR and Y ,v 1. 

X J 

The choice of a log scale for the axes makes the represen- 
tation simpler. A suitable unit for the plot is R. In 
J’igure 2 , OA = OB = OP = R. 


L L represents the plane IGHJ in Pigure 2. i 1/R 
X ■ y - 

on the line EP in this plane, and the limiit equation IT is 
■valid on this line. The domains of validity of the limit 
equations E and V are the portion above and belo^A? the line 
EE of the plane B„, that is EGHE and lEEJ respectively, 

excluding the line EE in both cases. 


As eq^uation (2.1) is symmetric in x and y, the order 

space plot is symmetric about the pilane I,. Consider 

X y 

the case when 1 cc L . The order of the inertia term is 
y X 

2 5 4- 

i that of the viscous term is '1 /Rl, . When 

the two terms are of the same orders Y 1= /l m;- E. which 

represents the plaiie EKLE due to the choice of log scales 

for the axes. This plane, to be called B^_, meets the plane 

J 

Ij„ nj along the line EE. Since 1 c. 1 , B is restricted 
X y ° y x’ y 

to only one side of the plane I . -..J L . limit eauation B 

X y y 

is valid on this plane except for the line EE. The dom.ain 
of validity of the limit equation E and V, are- regions 


above and below the plane B. such that 1 


that is. 


regions between the planes EGHE and EKLE and between the 
planes EIJE and EICLE respectively. A similar construction 
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is carried out when l^Ci, Ba.ch plane in the figure is 

X j 

infinite or semi-infinite. The construction of the figure 
is restricted tol 1i1 for clarity. 

The line joining 0 and Q in Figure 2 is in rhe plane 

L 1 and such that u = vu 1 . OQ meets the line 

X i y y 

of intersection PC of the plane L, 1 and the plane 

IK. 

B (EKLF) , at Q. Therefore the point Q represents the 

y 

classical hound ary- layer limit. Order coordinates of the 

__.i 1 

point Q are (1, R ^5 R ^). 

It is clearly seen from Figure 2 which shows the 

possible limits of NSEj that although there are infinitely 

many limit processes (each point of the figure corresponds 

to a limit process) the possible outcomes are finite in 

number. The domain of -validity of the limit equations E, V, 

B and B are planes wherea.s those for E , E , Y and Y are 

three-dimensional regions. Bach of these regions is bounded 

by two of the four planes EGKF, EKLF, BIJF, and Ei.ffiTF. For 

example, the domain of "'validity of Y^ is the region between 

the planes EIvMF and EIJF. Therefore limit equations E^, E^^, 

V and Y are each contained in the limit eQu.aticns corres- 
X y 

ponding to the planes vfnich bound their respective domains 
of validity. That is Y is contahned in Y and B and so on. 

X X 

The four planes mentioned above meet a.long the line EF 
which corresponds to RSE. The limit equations corresponding 
to these planes are therefore principal lim.it equations of 
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NSE. B and B. differ only in the definition of the 
X y 

coordinate axes. Thus we have only three principal limit- 
equations of hSE. They are the Euler equation, the Stokes 
equation and a third equation with which we associate the 
classical hound ary- laye r . However, this limit equation is 
more general. 

Ihe above construction of the order space helps in 

visualising various limits of HSE. A somewhat different 

construction of the order space can he one where we choose 

the order axes to he 1 ,1 and instead oft/ the order of 

X y 

the streatiF/ise velocity. This, however, represents a 
linear transformation in the order space. 

2.3 HiYSICAI FROBIEM 

The flow is considered to be confijied between two 
v;alls. A simple example of an internal separated flovj is 
the flow in a channel having a s5mimetric sudden expansion 
(Figure 1). Recirculating regions ore formed downstream of 
the expansion, and flow reattaches some distance downstream 
and finally evolves into a Poiseuille flow. An essential 
feature of this problem is that the flow is bounded in the 
transverse direction and, therefore, this length scale can 
be taken to be of order one. Two types of internal 
separated flows are identified later in Section 2.6. 

In external separated flo^-v, the interplay of the 



12 


recirculating flow and the enter irrotational flow requires 
that the shape of the hulihle or the hound ary of the 
irrotational flow he determined as a part of the prohlem. 

As a result the transverse length scale of the separated 
wake is difficult to ascertain. How'ever , the hound ar^y 
conditions on the two walls are sufficient to formulate the 
present r>T o hie m , 

2.4 LIMIT BQUAl'IONS WHEH L ret 1 

¥heri L .-v- 1, (i.e. the plane normal to L axis in 

y y 

Figure 2) possible limit equations with their domain of 

validity are shown in Figune 3. The limit equation K is 

valid at P. E and V are valid on the li' axis above and below 

P resrectively. B is valid on PA. whereas B on PD. Domains 

of validity of E., E„, and Y are regions OPA, OPDj APT 

j ^ y X 

and DPT respect iveljn Line OG in Figure 3 is such that 
V = - along the L^ axis u 1. There- 

fore, if the flow velocities are hounded for large Pl, the 
region of interest is that below AOG in Figure 3. Point G 
corresponds to the usual hound ary- layer w'hen the y axis is 
along the surface . 

For the physical problems considered the total flovj 
across any section remains constant and in its non-dimen- 
sional form is independent of R. Therefore , b'hv 1. From 
Figure 3 we see that there are then two principal limits. 
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They correspond to point 0 and point A. 

The limit equation corresponding to point 0(1, 1y1) 
is the Euler equation 

J ( y ) = 0 (2.3) 

The flov«i in this limit is inviscid . The above equation is 
elliptic in nature. Consequently , if the solution is 
disturbed at any point, its effect is felt at all other 
points of the region. 


The limit equation correspond ing to poiiit A (R ,1 ,1) 
in Eigure 3 is 


MJ 

ty 


ryps: 


'vi'' 1.-' 

'X iyyj' 



■ yj' IT 


(2.4) 


The length scale in the x -direction is of order R. The 
above equation is psnabolic so that there is a direction at 
every point, called the characteristic direction. The 
effect of a perturbation introduced at any point in the 
region is sensed along this local direction. 

2.5 DEVEIOHERT OR POISBUILIE EOT 


We first consider the development of Eoiseuille flow 
as a prelude to the internal separated flow problems, as 
the flow b.etYJeen the two par.allel vJalls develops ultimately 
into Eoiseuille f lo^-v . Wilson (1969) has considered the 
development of Eoiseuille flow using USE. A small 
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pei’fcux'ba.t ioi'i ana-lysis is caii^ied out; with the pelt urhat ion 
8,ssumed in the forni ^-^(y) exp(-Ax )5 wheie x and y aie the 
streamwise and transverse coordinates. An eigenvalue 
■problem is obtained for , X "being the eigenvalue. For 
la.rge E, two sequences of eigenvalues axe obtained. The 
members of the first sequence arc of order one, while those 
of the second are of order l/R, An eigenvalue of order one 
indicates that the perturbation is effective over a distance 
of order one. The eigenvalues of this clans are positive 
as v^ell as negative. The perturbations are, therefore, 
sensed in both the upstreaim and the downstresiin directions 
over a distance of order one. The streamwise extent over 
which the perturbations corresponding to the second sequence 
are felt is of order R. These eigenvalues are such that 
the perturbations decay in only the do"Vinstream direction. 

It then follows that if a small disturbance is 
introduced on the Poise'uille flow, it propogates both 
upstream as well as downstream. Eowter, the extent of 
upstream propagation is of order one and that downstreani is 
of order R. Also, it is the secorid mechanism which is 
dominant in the downstream propegration of any disturbance as 
its decay rate is considerably slower. 

The eigenvalue problem for the second type of eigen- 
function can be derived from the limit equation (2.4). 
Qualitatively the mechanism -of propagation of the 
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distuibaiice should therofore Tdg siniilc.1: to that of equa.tion 
(2.4). Since the equation b-s a, characteristic direction, 
the disturhence s are expected to propagate in only one 
d ire ct ion. 

It can then be argued on the basis of the nature and 
length scales of limit equations, and the analysis of Wilson 
that there are tvKj modes by vinich disturbances are propa- 
gated through the flovn If there is an upstreajn influence 
then this flovj is expected to be governed by the limit 
equation (2.3). Ihe second, mechanism is dominant for down- 
stream flow, and is described b3^ the limit equation (2.4). 

2.6 TWO TYIES OF ETTSffiUL SiriEiTIE FIOWS 

Two simple geometrical arrangements causing separation 
in internal flows are sudden contraction and sudden expan- 
sion (see Figujre 4). The separated flows in these configu- 
rations are qualitatively different. In the case of sudden 
contraction, the flow separates upstream of the geometrica.1 
configuration responsible for it. The streamwise length 
scale is, therefore, of order one , and the separation bubble 
appearing ahead of the contrantio.n will be of order one. 

The second type of separated flow is that downstream 
of a sudden expaiision. The dominant mode of propagation of 
the effect of the sudden expansion in the downstream 
direction is that corresponding to equation (2.4). The 
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streamwise leiigth scale is of order B. , thajt is a iyiDical 
lengths say, the length of the separated 1)111)116, hehaives as 
R. Viscous diffusion over a distance R is significant . 
Therefore, viscous forces are not negligible in the djanamics 
of the flow in the bubble. The 8.ppropriate limit equation 
is (2.4). Separated flov/s of this class are dealt, with in 
t he nr e se nt work. 


The above rather intuitive conclusion can also be 

reached by the following airgument . Vhen L 1 and 1 1, 

X y 

viscous terms are of higher order. In the 'near' limit, the 
flowj is governed by the Euler equation. Discontinuous 
surfaces ere allowed. One such surface is along the sepa- 
rating streamline. A vorticity layer is required for the 
vorticity to diffuse li the bubble. The thiclcness of this 
layer is of order (x/R) me therefore for any fixed x, 
it shrinks to zero in the limit. T'hc.n x R, the vorticity 
layer thickness is of the order of the transverse length 
scale of the charmel itself. The 'near' limit is, therefore, 
inapplicable when -x R. A 'far' limit given by R is 

require cl wh ic h le ad s t o e quat ion (2.4). 

2.7 EARLIER FRlSulTS 01 STPMI?riSE lENGTH SCALE 

Vi'e exam_ine the behaviour of a typical stream’vise 
distance for large Reynolds number in the earlier numerical 
calculations w;ith a, view to test the a.bove conclusion 
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regarding the nature of separated flow. 


H'orig and Macagno (1966) have calculated internal 
separated flow in a, channel having a syniinetric sudden ex- 
pansion. Their calculation is hased on JiSB and uses a 
f inite-d ifference scheme. Two methods have been used to 
study the prohlem. In the first, steadi^ state equations 
are solved, while the second appro a.ch is Dased on un stead;/ 
equations and the calculation is continued till a steady 
state is reached. The second approach is considered to have 
advantages over the first for large R. The ratio of expan- 
sion (downstream to upstream width) is two and the calcula- 
tion is performed for S (based on the average velocity and 
the half-width of the channel downstream of the expansion) 
upto about 160. Blow slightly upstream of the expansion is 
taken to be fully developed i.e. a Poiseuilla floxe. Two 


t;^ical streamwise lengths are shown as functions of R in 
Fig'ure 5. The;/ ane , the distance of the point of reattach- 
ment and the distance of the centre of the ’eddy' measured 
from the expansion. They cleanly show a line air trend for 
sufficiently lamge E. This problem is also considered by 
Morihara (l972). Calculations of KSE are performed upto 
R = 20, using a different numerical scheme. The length of 
the 'eddy' obtanned by iforihara is also shown in the figure. 
His results also show a linear trend. Flow in c gradually 
vndeiiing channel ha£ been calculated by Dorodnitsyn and 



shovjs the length of the 
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Me Her (1970). Their Pigur-e 6 
reversed flow region which increases 
linearly with R. 


approx imctel3' 


We thus conclude that the earlier calculations ha.sed 
on iTSE support the view tha.t the stressinwise length sc8,le is . 
of order R and it is thon expected that equation (2.4-) 
should govern the large Reynolds numher flow except possibly 
for some small sub-regions. 


2.8 EORMuMTIOE OE THE IROBIEM 


¥e shall consider the flow in a sudden expansion as a 
prototjrpe lor our discussion, a.lthoagh the conclusions exe 
expected to be more generaR and applicable to other similar 
problems such as shown in figure 1. 

The flow downstream of the expaxision is die'ided into 
two sub-regions corresponding to H-'V- 1 and R. The 

limit equation for the ’far’ region (l^-'^nR) is (2.4). As 
a result of the limiting process, three derivatives in x are 
lost. Equation (2.4) can therefore satisf^a onlj^ one condi- 
tion in X vvhich is obtained on ma,tching with the solution for 
the ’near’ region 1). The ’near' region solution is 

somewhat difiicult to construct. As the limit equation for 
this region is elliptic in nature , it would further require 
matching with the far upstream flow. To obtain the initial 
condition for equation (2.4) in the absence of matching, we 
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maice a \^jorl<:irxg hypothesis’ -^^hich is to be tested subsequently, 
thc^t the initial velocity profile is the same as the far 
upstream flops'. 


We no?J eliminate R in equation (2.4) by scaling the 
X coordinate b3^ R. It then becomes 


'I' 


! y 




'4^' ii-'* 

nx lyyy 


7JJJ 


(2.5) 


Unless othervv'ise stated , x Vvould henceforth implj^ the 
contracted strsam\?.‘isG coord inete. In this work E is taken 
to be based on the av/erage velocity end the half -width of 
the channel downstream of the expansion. 


The initial condition required for the integra,tion of 
equation (2.5) is given as 

^-4x = 0, y) = • (2.6) 

This condition is equivalent to prescribing the strearm^ise 
velocity at the initial section. The no-slip condition at 
the wa,lls requires 

r'= ± n ’iv =0 y = ± ^ • (2-7) 

■ J 

In the coordinate sj^^stem x (contracted) and y, the 
'near' region shrinks to a plane, x = 0. Equation (2.5) can 
then be regarded as the field equation valid for x> 0. ft 
X = 0 an initial condition (2.6) is prescribed. Equation 
(2.5) along with conditions (2.6) and (2.7) constitute 
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the ma.themat ic al probleni aiicl subsequeiitly vie test ii they 
form a. sufficient basis to calculate the flow. 

Equation (2,5) can he integrated once with respect to 
y. The constant of integration is the pressure gradient. 
However, as we wnll see, it is preferable to work wutli the 
higher order equation (2.5) than v-? it h its integrated version. 

2 . 9 EHRTHEEi COMiSENl S 


The equation obtained for the calculation of flow 

development is similar to Prandtl’s bound aj!:3’'-la3rer equation 

BLE. However, it should be noted that the viscous region 

has a. thickness of order one whereas the thiclaiess of the 

_ 1 

classical boimd a.r3""laiyer is of order . Also, instead of 
pressure being impressed by the external flow, it has to be 
obtained as a pent of the calculation here. Eow;ever3 there 


are four boundan37' conditions in y to be satisfied. 


.ne 


confinement of the flow can thus be said to generate its own 
pressure field. Strong interaction of the separated region 
with the mainflov/ wihicb is a feature of the external 
separated flow? is absent in the present internal separated 
flow . 


It ma3^ appear surprising that the equation which w/e 
intend to use for flo’WB having recirciilait ing regions is of 
the bouuid ar3''-lei3?'ei type. BliE are knownn to pjossess a singu- 
larit3- at the point of sere skin friction (G-old ste in , - 1 9.48)) , 
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aiid it is "believed that the continuation of the calculation 
beyond this point is generally not possible (Brovm & 
Stewajctson, 1969). Flovj reversal occurs beyond the sepa- 
ration point. Since the argument used in deriving BIE 
does not require the absence of reversed f lov; , ElE permit 
reversed flow a.& long as the bound layer thickness is 
of order R~^. Briley (1971) obtained the BS solution for 
the flow in a separation bubble by means of a finite- 
difference scheme, Y\7hich with increasing R appears to 
approach the corresponding solution of BLE. 

Consider Figure 6 (taken from Brown 1 Stewartson, 
1969). I'he conditions on an initial section RB, the wall 
ilF and the edge of the boundary-layer BB are given. Since 
the flov; in GEF is reversed, the convective propagation of 
a disturbance in the sense of decreasing x is possible. 

In such circumstances, the pamabolic equation is said to 
be of the mixed type. There seems to be an intimate 
connection bet?^een the theory of paraibolic equations of the 
mixed type and separated f lov; . Apparent 13 ^, the requirements 
oil the boundary condition are severe for such equations. 

The theory of such equations has not yet been sufficientlj- 
d eve loped to give conditions which the given data must 
satisfy in order that the problem is well posed and a 
solution exists. In the absence of any formal theory, 
qualitative arguments are advanced for the approrr isite 
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"boundary condi'fcions. It may appear "bnat in addition bo ths 
already specified conditions on IB, BD and AB, a condition 
on EB (Bigure 6) is required in order to determine the 
solution over ABDB. Hovt/ever, if the toujidary conditions 
7 »ere to "be completely ar"bitrary, flow upto point G will "be 
affected by changing the condition over GF and EB, and 
possibly the separation point G itself will shift; but the 
flow in ABCG uniquely determines the point G. Perhaps some 
condition on EB is needed in general to determine the eigen- 
solutions if any. The above argument also suggests that the 
boundar-y condition to be prescribed should be in a certain 
class for the solution to be regular. The regularity of the 
solution imposes certain restrictions on the boundary/ condi- 
tion, even when the region of interest is A3CG. when the 
given condition is otherwise, the flow manifests a singu- 
lar ity. 

Based on the above discussion, the appearance of a 
singularity can be viewed from twn angles. If we are 
looking for 'thin' separation (the contention here is that 
it is regular) , the singularity represents the incompatible 
nature of the supplied data. Catherall and Mangier (1966) 
designed a thoughtful scheme in which the supplied data had 
partial freedom to adjust in order that the solution be 
continued in a regular \Na.y at each step. On the other hand , 
the singularity -would mean the failure of the bound ary-layer 
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limit as an approximation to t^e actual jjhysical phenomenon 
('thick’ separation). 

There are two distinct issues in the calculation of 
separated flow, namely, the equation and the hound ary 
conditions required for the formulation of the mathematical 
problem, and a. method of calculation for the formulated 
mathematical problem.. Equation (2.5) is considered to be 
appropriate for the present class of problems. The 
boundary conditions prescribed on the walls (2.7) seem to 
be of the simplest type. li any difficulty arises, it can 
be attributed to the initial condition (2.6). Yfe shall 
regard the problem as well posed and assume that a solution 
exists. 


The usual methods of integration of par’abolic equa- 
tion are not expected to be applicable in the presence of 
reversed flow. An iterative technique is used by Klemp and 
Icrivos ( 1972 ) in which the shape of the reversed flow 
region is assumed , and then the solution in the forward 
direction and bacto'Jard direction are obtained separately by 
a marching technique suitable for parabolic equations. The 
shape of the reversed flow region is iterated to match the 
flow quantities across it. In solving a parabolic equation 
of mixed type, Bai and Euerti (1969) have used a numerical 
scheme Vi/hich is normally applied to elliptic equations. 
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Integral relations methods or methods of weighted 
residuals are well knovm for solving partial differential 
equations. Their convergence can he shown under certain 
circumstances. The method of cad-culaticn used in the 
present work to solve equations (2.5 - 7) belongs to this 
class and is developed in the next chapter. Hov^ever , the 
approximating and weighting functions are eigenfunctions 
of a related problem, unlike the usual functions which are 
often selected for their analytical properties and 
computational Simplicity. 



CHAPTEE 5 

mmoD OP caeculatiow 

3.1 BITRODUCTION 

This chapter deals ■with the formulation of a method 
to calculate the flow de'velopment heti^ijeen two parallel walls, 
which can be applied to internal separated flov»;s. The 
governing equation and the required boundary conditions are 
given in Section 2.8. 

Since the flow evolves ultimately into Poiseuille 
flow, equation (2.5) can be linearised about this flow for 
sufficiently large downstream distances. It then leads to 
an eigenvalue problem and the perturbations decay exponen- 
tially with increasing streamwise distance. In the region 
where the flow differs markedly from Poiseuille f lov? , the 
linear approximation does not hold , and nonlinear terms must 
be retained. Ko'wever, if the nonlinear interaction terms 
between the modes deca^'’ rapidly, an expansion in the eigen- 
functions could lead to an efficient computing proced'ure , 

5.2 LUEAR PEOBIEM 

As x-^oo , the channel flow develops into Poiseuille 


flow given by 
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Let the flow he expressed as a sum of loiseuille 
flow and a perturbation , given by 


•p i1 


(3.2) 


Substituting (3.2) in equation (2.5), we get 
kU 

■lys'yy iP nyyx 


+ ''■id H^iy '^iyyx ~ tix '^yyv'^^*^^ 

ip !lyyx ^ p 'lx ' 


Lash indicates differentiation with respect to the argument. 
Neglecting the quadratic terms in ^ and substituting for 
from (3.1), we get the linearised equation 

'i'liTyy ~ - y ) 'Viypc " ^ '"Hx ” ^ * (5-4) 

The form of the above equation suggests a solution of the 
form 

= <P(7) , (3.5) 

for some X . Substituting, (3.5) in (3.4) we get 

+ X [ 3 / 2(1 - y^) o + 3^] = 0 , (3.6a) 

where X is an eigenvalue. The wall conditions require that 

4^ = o' = 0 at y = + 1 . (3.6b) 

Linearization aLout the Poiseuille flow thus leads to an 
eigenvalue problem EVP, and the perturbations decay expo- 
nentially in the positive x direction, if X is positive. 
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5.3 'IHS EIGEKTiiUE iEOELEM 

The form of (3.6) suggests that its solutions csn be 
separated into odd and even solutions. An eigenfunction is 
called odd or even depending on the antisymmetry or the 
symmetry of the associated stream'wise velocity/ component. 

EVP (3.6) has been studied earlier by Gillis and 
Brandt (1964) and by Wilson (1969). Gillis and Brandt have 
made a direct numerical computation of even eigenfunctions. 
Wilson has obtained equation (3.6) as the equation governing 
the leading term in a large Reynolds number expansion. 

Equation (3.6) can be integrated once to give 

h + X '|3/2(1 - y^) O + 3yX|= Constant . (3.7) 

The constant of integration in (3.7), u’ben multiplied by 
e" gives the associated , pressure gradient. Evaluating 
the left side of (3.7) at y = 0, the constant of integra- 
tion turns out to be zero for odd e igenfi.inct ions. Therefore 
the pressure field associated with an odd eigenfunction is 
uniform. 

It has been shown by Wilson that the odd eigenvalues 
are real. It has not been possible to prove a similar 
result for even eigenvalues. Hov/ever all the eigenvalues 
calculated by using numerical methods are real (Wilson, 1969). 
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The real eigenvalues of (3.6) can be shown to be positive 


by multiplying (3.6) by 4^ and 
leads to 

'■1.2 . r -' 

dy - 3/2 A j i (1 - 
il 1.-1 


integrating by parts, which 


2 '2 
y ) <^ <3y 


1 

j dy 

-1 


= 0 . ( 3 . 8 ) 


Now consider 

Z' , 2 1 ^ 

X = i (1 “ y^) 4^' dy/ (4:^^dy 
~1 ' 

Stationary value of leads to Legendre's equation. The 
minimum value of , therefore, is two. Hence the term in 
the bracket of (3.8) is positive. Therefore ,X is positive. 


t'ilson has reported that manhine calculations have 
shown no real eigenvalues in the range -200 <. X < 0. 
Although it has not been definitively established that all 
the eigenvalues are real, it is conjectuT'ed that they are 
real and therefore positive. 

Appendix 1 contains an outline of a scheme for 
calculation of eigenvalues and eigeniunccions . fourth order 
Runge-Kutta-Gills method has been used to integrate the 
differential equation (3.6). Eigenfunctions of higher order 
oscillate more rapidly and consequently have larger deriva- 
tives, which require a suitable method of integration. 
Estimates of A for higher modes have been obtained in 
Section 3.3.1. Table 1 contains the first few eigenvalues. 
Thej’- are numbered in increasing order and are alternatively 
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odd and even. The lowest eigenvalue is odd in nature, and 
hence the velocity perturbation which persists the longest 
is anti symmetrical. 


m 

1 

2 

3 

4 

5 


m 

14.4535 

18.8146 

48.8857 

57.5463 

104.572 


TABIE 1 EIGEI'TVALUES OE EQUATION (3.6) 


3.3.1 ESTEvIATES EOR LARGE Ox 

Analysis of equation (3.6) for large X is soniewhat 

similar to that of the Orr-Sommerfeld equation. There are 

four independent solutions, two of which correspond to 

length scale one, and are denoted by <f- and h . The 

^ "'E 

differential equation governing ■^op. and , in the 

■'A -R 

limit cc , is 

(1 - y^) cto + 2cr. = 0 . (3.9) 

One of the solutions, is a regular function. However 

is singular at the w;all (log singularity) and therefore 
an appropriate ’wall solution’ is required. 
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The other two solutions are rap id ly varying . These 
solutions are obtained using the V/KBJ method . he denote 
these solutions by and 4'V. • sxe obtained by 

writing 


^'g,d ^ • 

We assume a series solution of g in powers of X ^ 


( 3 . 10 ) 


g 





+ g. 




(3.11) 


The series solution so obtained is 

X. 

B = ± ^ ^ (1-y^)^ + 5/2 y(l-y^) + ... , (3.12) 

ViJhich is valid only if 


- yhh>i 


(3.13) 


Substituting (3.12) in (3.1C)j we get 


2^ ^ ^ ^ Cos il^■3/o'il^2 II A \ , £-<■; ’• ! -1 / > 


%,L ^^-y ^ Si5 i 2 ( 3 / 2 A )2 (y ^ gin y)/.(3.14) 


Condition (3.13) is violated near the walls and therefore 
the above solutions are invalid there. The proper stretch- 
ing for the solutions near the wall can be obtained by 
considering the limit of (3.14) when 3.' approaches the wall, 
let 


s = 1 + y 


(3,15) 


Then as s 0, the limit of (3.I4) is 
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(-^/2 + 4/2/3. s^/^)J. (3.16) 

The proper stretching required near the wall, therefore, is 
-1 /3 

X . This can also he obtained by considering equation 
(3.6) or (3.13). The 'near wall' solutions are obtained in 
terms of Airy's function and matched vvith the WiCBJ solutions, 
which gives 

^.^rvin^ . (3.17) 

Detailed calculation leads to 

Odd j X^--'32/3 (n + 1/12)^ , (3.18a) 

Even j >.^--'32/3 (n + 1/4)^ . (3.18b) 

It should be noted that the leading term of eigen- 
values are real and positive, and the eigenvalues are 
countably infinite. Dor a given n, the odd eigenvalue is 
smaller tha,n the even. If the eigenvalues are arranged in 
the increasing order, they ar-e alternately add and even. 


3.3.2 ORTEOGOIIALITY EILAT IONS 


The present eigenvalue problem is of the form 


uA. 4>+ A jV\ X = 0 , 


■X' o = 0 


(3.19) 

( 3 . 20 ) 


. . . 

and are the linear operator 


where 
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^ = dV(3y^ , vH.-- = 3/2(1 - y^) + 3 s (3.21a) 


and ^ represents the hound ary conditions 


<f'' = 0 ; 4 = <4 = 0at3r=:+l 


(3.21h) 


The eigenvalue lorohlem encountered here is not of the usual 
type. It can he considered to he of a more general nature 
and reduces to an ordinary eigenvalue problem, v'henjX/j is 
one. vJe shall assume that the eigenvalues are real and 
therefore positive and denumerahly infinite, how a scalar 
product is defined hy 


_1 

( cb, Q ) = i oGdy. 

-1 


( 3 . 22 ) 


Corresponding to and .AIj > define the adjoint operators 
hy the relation 


4 = dP^(<p,ii )/dy , 

rp>; = dP 2 ( 4,'r )/dy , 


where th and 4 are any two suitably continuous functions 
and P^ ( 4- , o ) and p 2 ( 4o , B ) are bilinear concomitants. 

By suitably choosing the boundary conditions so the.t and 
p 2 vanish at the boundaries, we state the adjoint problem as 

= 0 . (3.23) 

Boundary conditions for the adjoint problem turns out to be 
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the same as (3.20) and the adjoint operators are 

= ^ = dVdy"^ , 

Jd.-* = 3/2(1 - y^) d^/dy^ - 6y d/dy , (3.24) 


If and are e igenfimct ions of (3.19) and 

(3.23) respectively'" and the corresponding eigenvalues are 

A and /\ 5 it can he shown that 

m n 


( -A 


m 


dN ) (/'::.u , 
n ^ ^ ■ m 


e,j = 0 


■Ihen m / n, it then follows 



i-4 ) = ( , 7" P, ' 

'-n m ^ n^ 


0 


or 0^ dy = 0 , m / n (3.25) 

lie can define a new scalar product ?Jith respect to the 
operator X. j "by "the relation 


-+ 5 a i 
,toj. 


= ( tC hc>, y ) 


(3.26a) 


where dt and '0 satisfy (3.2 0) and ore" suff iciently "d ifferen- 
tia,hle. Iherefores 


Equation (3.25) then becomes 

, P) i = 0 , m / n 
L m ’ '^ni 

Eigenfunctions of (3-19) and of its adjoint 


( 3 . 26 b) 


(3.27) 

problem (3.23), 
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therefore 5 form a h fort hog onal set, where the orthogonality 
relation is given hy (3.27). 

3.3.3 EIGEIOWCTIOIT EXPih^SIOl 

We Yoill require the expansion of functions satisfying 
suitable conditions into eigenfunctions of the equation 
(3.19). We shall assume that the required eigenfiunction 
expansion converges in some sense. 



The perturhat ion stream function is small 

sufficiently downstream, so that the nonlinear interaction 
of the modes in (3.3) can he neglected. The solution of the 
linear problem (3.4) can then be written as 


(x 

rn 


) 


(3.30) 


The constant A in the expansion of h . can be found using 
the relation (3.29), 
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A 


m 






9 


where 


(3.31) 


. - = '4-' at X = X 

10 1 o 

3,4 mmOD OF CAICinATION 


(3.32) 


let in equation (3.2) he written as 

= X a^,(2^) 4^01 (y) > (3.33) 

m=1 

where the coefficients depend on x, and 4“^^ are eigen- 
functions of equation (3.6). In an H-term calculation 
summation index in (3.33) runs from one to H. The range of 
summation of indices will not he henceforth mentioned as it 
would he clear from the context. Substituting (3.33) in 
equation (3.3) followed hy multiplication hy , and 

integrating, we get 


! 


3. + 


:\ 

m 


a 




>'.T;c 

fr q mpi 


a a„ 
P 1 


where 


(3.34) 


^mpq "^m ^impq^^Sm ’ (5^359.) 

^Impq “ ^ "'m ’ '"^p ^q ~ ‘^'p '^q ^ ’ (5.351) 


r 


"2m 


4 ft 
_ ■ m ’ mj 


(3.35c) 


The differential equation (3.34) gowerning the coefficients 
of the eigenfunction in the expansion (3.33) is a set of 
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first order ordinary, quasi-line ar , autonomous differential 

f 

equations. The first term on the left in equation (5.34) 
represents the linearized part of the inertia term vvhereas 
the transverse viscous diffusion is represented hy the 
second term. If the nonlinear interaction of the modes 
given by the term on the right side of (3.34) is negligible, 
a linear equation is obtained which has the solution (3.30). 
Equation (3.34) can also be ^’’■jritten in the form 


/ ' ( h “VC a ) a = - X a 
mn mpq a ^ p mm 

p - - 


(3.36) 


The inertia terms have been collected on the left side ,in 
the above form. is the Kronecker delta. 


3.4.1 IhITlAl CONDITION 


The flow calculation is now reduced to solving 
equation (3.36) with appropriate initial conditions. 
Initial values of a^^ can be obtained in a. similar via:y from 
(3.33). At X = 0 

tio = V • 

Therefore 


a 


mO 


a^_(0) - ( 


. IV 

’m 


) / I 


2m 


(3.37) 


If prescribed in equation (2.6) is piecewise t?^ ice 

differentiable, so is V^q* 


Equation (3-37) then leads to 
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^mO = 


i ” 


e. ) / X 


m 


2m 


3.4.2 LDTEAR S'UB-EAi^GB 


(3.38) 


As X goes to infinity, coefficients go to zero. 
However, the computation need not he carried heyond a certain 
stage when the coefficients are small, in which case the 
solution of the linea.r problem can he utilised. Then 

a = a (x ) exp | - /■ (x - x ) ! , ( 3 . 39 ) 

m m o m Oj 

where x is the streamrjise coordinate heyond vvhich the non- 
o 

linear term in ( 3 . 34 ) are negligible. 

3.4.3 SmASIRIC PROBIEM 

When dealing with a S 3 nrimctric problem, only even eigen- 
functions need to he considered. Index m in (3.33) then 
runs over the set of U even e igenf met ions . The range of 
integration with respect to y in any integral such as ( 3 . 22 ) 
can therefore he restricted to the upper half of the 
channel, 0 3 ' Ar. 1 

3 . 4.4 Om -TERM SO lUT ION 

When onlj?- one term is used in the expansion (3.33) the 
resulting differential equation (3.36) can he integrated to 
give 

ln(a^/a^Q) - 0^^^(a^ - b^q) = - '-^x , (3.40) 

where a. p = a. at x = 0. As C = 0 for a 1-term 
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asymmetric calculation, vi/e oltain 

- A.x 

^ ® ^ 

3.5 OH EQUAIIOH (3.36) 

Equation (3.36) obtained for the calculation of flow 
development can be ’/written in the form 

C a = - /'i a s ( 3 . 41 ) 

where a is a vector whose elements are a^ , a2 > .... , 41 is 
the eigenvalue matrix, and C may be termed as an inertia 
matrix. C and A ajce H x H matrices and a, is a H vector in 
a calculation based on IT terms. 

Equation (3.41) is a set of algebraic equations for 
a , which can be solved using any standard method, A 
typical integration scheme would then be as follows. At a 
given X, the derivative a is obtained from the given a. 
Integration to the next point in x is then carried out using, 
say, t he Eule r me t hod . 

The derivative a’ can be calculated provided C can 
be inverted, which is possible, only if the determinant of 
C is non-zero. The integration procedure, therefore, 
breaks down for the chosen value of H in the event the 
determinant of C becomes zero at any x. The singularity of 
the resulting equation (3.41), which excises ’when the 
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determinant of C is zero, is a general feature of the 
equations of the alove type. Such a difficulty can occur 
when a moment method or an integral relations method is used 
to reduce partial differential equations to ordinary differ- 
ential equations similar to (3.41). We take up a discussion 
of the pha.se space to examine fur-ther the nature of the 
singularity. 

3.5.1 PHASE SPACE 

An N-term solution can be represented by .an integral 
cur've in an H-d imensional i^hase space, whose coordinates 

are a^ , a 2 > .a^-. The equa.tion governing such integre'L 

curves is 

da.,/l).| = da2/D2 = = daj^/f.- , (3.42) 

where D is the deteraiinant of the matrix olteined 'bv 
mi 

replaning the m^th column of the matrix C oy the vector on 

the right side of equation (3.41). As equation (3.35) is 

autonomous, E and D, the determinant of C, are functions of 
m 

a^. P'urthermore , they ar'e poljrnomials of degree U. 

The origin of the phase space corresponds to the 
fully developed flow and therefore a ' meaningful' solution- 
should approach this point. The given initial condition 
determines the integral curve in the phase sp.ace . The 
evolution of a with respect to x which can he reg:arded as 



40 


a parameter of the integral curve obtained by 

solving one of the equations 

The condition 1=0 would normally represent a sur- 
face in the phase space . If an integral curve crosses the 
surface 1=0 (Pigure 1 ) , the direction of increasing x 
appears either 'sink' like ur 'source' like and is ph 3 /'si- 
cally unacceptable for describing a real flow. However, aio 
integral curve could cross the surface 1 = 0 at certain 
singular points of the equation (3.42), if there are any. 

The phase space of symmetric calculation for N = 2 shows 
such a po£ sibilitj)-, and is considered here as an illustra- 
t ion. 

3.5-2 PHASE PLAilE' 

The axes of the phase plane are a^ and ag , and the 
differential equation governing the integral curves is 

da^/dag = * (3.44) 

Explicit forms of 1^, I 2 , and 1 are given in Appendix 2. 
Conic sections that are associated with the quadratic 
expressions, 1^ and I 2 , are ellipses and with 1 is a 
hj^perbola. Eigure 8 shows the ellipses 1^ = 0 and I 2 = 0. 
Evidently, when an integral curve crosses the ellipse 
1^ = 0 or I 2 = 0, except at their intersection, it is 
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normal to the or the a 2 axis which is shown hy "vertical 
and horizontal arrows in Figure 8. 

Both = 0 and = 0 pass through the origin 0 and 

are normal to a^ and a„ axes at 0. Iherefore , 0 is a 
12 

singular point of (3.44) j and it can easily oe shovm that 0 
is a sta'Dle tvjo-tangent node with a^ axis as the tangent. 

Figure 8 shows another point of intersection of 

D = 0 and 3 =0, which is denoted oy P. Equation (3.44) 

1 ^ 

therefore has another singular point P which is shown to he 
a saddle point in Apoendix 3. 

The circve D = 0 is also shown in Figure 8. The curve 
is seen to pass through P. This can also he shown from the 
differential equation. If D.| =0 and D 2 = 0 have a point 
of intersection other than 0, at which D 4 ^ Cq then from 

II 

( 3 . 43 ) a^ and a^ are zero. The left side of (3.36) is 
therefore zero, whereas the right side is not. Therefore, 
we conclude that D = 0 at P. 

Given an initial condition, the integral curve in the 
phase plane can he traced. A typical integral curve is 
shown in Figure 8. If the integral cunve crosses I) = 0, 
except at the second singular point P, the solution is 
physically unacceptable. However at P, as we shall see 
below, it is possible for an integral curve to cross from 
one side of the curve B 


= 0 to the other. 
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¥e consider a schematic of a saddle point singularity 
in jpigure 9, where hollow arrows show the direction of x 
increasing if D were positive everywhere, low let us draw 
a curve D = 0 passing through the saddle point, let 1 0 

above the curve D = 0, and 1 <0 "below it. Ihe direction of 
X increasing is now shown "b;/ full arrows. He see that there 
is only one wa3' of crossing from a given side of the curve 
D = 0 to the other that is phj^sicallj^ acceptable. Figure 9 
also shows the 'sounce' or 'sink' like "behaviour of the 
integral curves a,t points other than P of the curve 1=0. 

A point in the phase space which can "be joined to the 
origin an integral curve which is physically,' acceptable 
will be called accessible. The set of accessible points 
will be called the accessible domain. 

It is therefore seen that all points of the phase plane 
are not accessible. If the given initial condition for 
equation (3.36) is such that fne initial point is accessible s 
then no difficulty is encountered in the integration of 
equation (3.36). 


3.5.3 CONCLUIIHG BEMARKS 

It is seen in the previous section that some of the 
points of the phase plane are not accessible due to the 
singularity of the differential equation. The occurttmce of 
I) = 0 and the presence of the singular points are the 
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features associated vv'ith the equation. They can. therefore 
occur, in general, with any number of terms in a calcula- 
tion. lor example, in the case of 1~term sj-'mrnetric 
calculations, the accessible domain is a.. 1/Ch In. - 2.5. 

3.6 THE DTTEGRATIOH PSOCSDURS 

An integration proced'ore which proceeds in the phase 
space, is adopted with a ‘view to have a simple control over 
the accuracy and the speed. Usual methods of integra,tion 

C 

with X as the independent variable will be applicable when- 
ever the present method is successful. At each step of the 
integration of the differential equation (3.36), the 
derivatives 


8 = E , m = 1 , 2 , 

m m 


(3.45 




are determined, let n be such that 

, m = 1,2, , U . (3.46) 

Then a is chosen to form the independent variable for the 
n ■ 

next step of integration. Equation (3.36) can then be 
written as 

" \^®n ’ (3.47a) 

dx/da^^ = • (3.471) 

A simple pred ictor -corrector method is designed, to 
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solve (3.47a). It maintains a uniform level of accuracy 
throughout the calculation. This is effected hy incorpo- 
rating a method of reducing the step length in a^^^j to 
achieve the prescribed level of accuracy" when the deriva- 
tives are large, and on the other hand, increase the step 
length to speed up the computation whenever possible . The 
method is self starting. It is described in Appendix 4 
which also includes a brief description of the computer 
progra.m. 


As discussed in the previous section, equation (3.36) 
is singular: at 1 = 0. Therefore, if during the integraf;ion, 
D vanishes at any stage, the integration is terminated. 

The above integration proced'ure is used upto a point 
where the varia.bles a^ are small enough for the solution of 
the linear problem to be applicable. If ia^j C'^ifor all m, 
where Cm is a prescribed small number, the integration is 
terminated, and equation (3.3S) is used to extend the 
calc ula.t ion upto downstream infinity. 

3.7 RESUITS OF CAICUIATIOIJ 

The output of the calculation is primarily a^^(x). 

One can then determine several flow variables. For instance, 
centreline velocity u ; pressure gradient dp/dx and wall- 
vorticity o can be determined by 



45 


= ''•5 + , (3.48a) 

■dp/dx = -3 + (1) » (3.481D) 

^ '--I,. = ± 5 - .La^4''j^(+ 1) . (3.48c) 

The values of at the upper and the lower Y<alls are denoted 
hy suffix + and - respective!;'-. Stream f-unction is deter- 
mined using equations (3.2) and (3.33). 

Re attachment point on the wall is determined hy the 
location ?ihere weJ.l-vort icity changes sign. In the case 
•where reattachment is on the centreline , the chaiigs of sign 
of u^ determines its location. The centre of the eddy is 
determined by finding the location of the extremum of 
stream function within the flow. This operation also yields 
the vai,lue of the recirculation. 

In addition, tv?o quantities which are of general 
interest are considered below. 

3.7.1 f S'VE 10 imm lENGTH 

to 

Development length is ana-logous^-ahe entry length and 
iT. can be defined as the streamwise distance after which the 
centreline velocity is within a given fraction of the 


ultimate value. 

for example , 




^99 “ i 

%(xQg) - 1-5i 

/1 .5 

II 

o 

b 

(5.49a) 

^999 ° ' 

^0^^999^ - ^*51 

/1 .5 

= 0.001 . 

(3,49b) 
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That is, Xq (3 is the value of the streamwise coordinate at 

J J 

which centreline velocity u differs from its ultimate value 

0 

1 .5 hy '\'l- • Similarly corresponds to a difference of 

0 . 1 %. 

If X is a streamwise location after which the 
0 

solution (3.39) is applicable, then the development length 
can he found from 

i ■ V a (x ) exr ■- Tv (x - x ): l-'(O): = 0.01 x 1 .5 • 

The contr ihution of the modes higher than the first even 
mode, can usually he neglected, as they decay more rapidly. 
We then have 

>^99 = (0-01 X 1.5/|ag^(x^)j) ,(3.50a) 

where the subscript el denotes the first even mode. 

Similar ly 

^QOO (0.001 X 1 .5/ ;a„., (x )i ) .(3.50h) 

J J J 1 t C 1 O ; 

3.7.2 PEESSlll:E EI^COVERY COEI'FIC TEST 

The non-d imensional pressure graidient is --5 in 
Poiseuille flowj. We define p^.(x) as 

dp^ydx = dp/dx + 3 . (3.51) 

We taice p( O) = Pj,(0) = 0. The pressure recovery coefficient 
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is then defined "by 

) • (3.52) 

The asymptote to the curve p(x) is therefore given "by 

p = -3x + . (3.53) 

If the solution (3.39) is applicable beyond x = x , ^ 

o "r 

is gi-ven by 

S? t 

^ ^r1 ^ -'''m ’ (3.54) 

Y;he2:e 

Pr^^o^ = P^'^o^ + ^""o • (5.55) 



CHAPIER 4 

EKl'RY -PLOW DI A CHAKLEL s A TEST CASE 


4.1 lETLODUCTION 

We consider the entrj/ flow in a straight channel as 
a: first application of the calculation method developed in 
the previous chapter. 

Although the problem is technically important, the 
attention it has attracted is considerably more than what 
might he expected . Its simple geometry provides a suitable 
testing ground for theories and methods. The present 
application is aimed at establishing the method by comparing 
its results with those obtained ea,rlier by other methods. 

The comparison provides an indication not only of the 
accuracy of the method of integration but also the applica- 
bility of the limit equation and the initial condition 
hjj^pot he s is . 

There are a large number of reported investigations 
on this prob3.em. We, however, consider mainly those which 
are required for the comparison of the results. The first 
anal^'tical work based on suitable approximations is due to 
Cchlichting (1934> 1968). The flow near the inlet is 
assumed to consist of the wall boundary-layers together with 
an inviscid uniform central region which accelerates to 
satisfy the continuity requirement. A 


series method is 
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used to construct the ho unda.ry- layer solution in the initial 
portion. The solution for large downstream distances is 
obtained hy perturbing the Poiseuille flov;. The two 
solutions are patched at an intermediate location. This 
solution has been later improved by including additional 
terms in the perturbation series by Collins & Schowalter 
(1962, 1963) and by Roidt & Cess (1962). Bodoia & Osterle 
(1961) have obtained the solution of BIE using a finite- 
difference scheme. However, a large truncation error is 
unavoidable, further, the scheme is unstable in the 
presence of reversed flow. 

Hang & longwell (1964) have considered two types of 
entry conditions in their solution of complete USE obtained 
by a f inite-d iference method. The flow is assumed to lae 
uniform at the channel entry in the first case and a.t the 
upstream infinity in the second . The solution for the first 
case is also obtained by Gillis &, Brandt *^1964) and by 
Morihara (l972). in off-centre maximUftcf streanrAuse velo- 
city occurs after a short distance doYmstream of the inlet 
section, and it continues till it is eroded by the thicken- 
ing of the boundary-layers. This apparent 13^ paradoxical 
situation has been clarified in an important paper b’/- 
Van Dyke (1970), who also showed the singular perturbation 
nature of the entry flois problem. The problems considered 
by Schlichting (l934) and by Bodoia & Osterle (196I) are 
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essentially one-term dovaistream expansions. The first type 
of entry condition appears restricted and introduces frac- 
tional povjers of Reynolds number iii the expansion, 

4.2 THE IHTIAL COK'D ITICH 

The channel consists of tv;o plane parallel plates 
located af y = + 1, x >0 (Figure 10). As the flovj is 
symmetrical in y, only the upper half of the channel is 
considered. 

The flo\v calculation as formulated in Chapter 3 is 
reduced to sol-ving equation (3.36) v;ith an appropria^te 
initial condition. At the inlet, flovj is taken to be 

u = 1 , 0 4^y4l , ( 4 . 1 ) 

= 0 , y = 1 . 



The number of terms H to be used in a calculation is 
to be chosen in the first place. The ability to predict 
to a good approximation with a few terms is a great advantage. 
We intend to use rather small number of terms in our 
calculation. A 1-term solution is given by (3.40). 


51 


T'wo , three and five terms have heen used in the calculation 
whose results v^e discuss in the next section. Galculation 
with more terms is, in principle, straightforward. 


4.3 SESlJliTS 

She inlet velocity lorofile is approximated as shown 
in Figure 10. The given uniform velocity profile is approxi- 
mated hy functions which vanish at the . The region of 

significant error tends to get localised near the wall with 
increasing number of terms. This evidence can "be regarded 
as a coniputat ional indication of the convergence of the 
approximating functions, although the nomunif ormity in con- 
vergence neau the wall is to he expected. A similar situa- 
tion arises when a unit pulse is approximated by trigo- 
nometric functions. To balance the reduced mass flovj near 
the wall, velocity elsewhere has to Ice greafer than unity. 
Table 2 gives two indicators of the approximation, namely 
the centreline velocity u at x = 0 and "the associated 

CO 

momentum flux . The momentum flux is related to the inte- 
gral of the square of the error in the streamvise velocity, 

1 1 

Hu - 1)^ dy = I u^dy - 1 . 

o’ 0 

It is seen from Table 2 that both the quantities approach 
the given value v;ith increasing number of terms . 

I.i ... 

. „ 4 51152 ^ 

Ate. No, 
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N 

2 

3 

5 

Given 

u 

CO 

1 .206 

1 .122 

1 . 090 

1.0 

M 

0 

1 .088 

1 .071 

1 .051 

1 .0 


TABLE 2 REPEESEBTATIOK OF THE' ILITIAL rROFIIE 


Figiixe 11 sho'ws the calculated centreline velocity. 

As mentioned shove, u is initially different from unity. 

L/ 

A convergent trend is clearly’- seen from the figure. In 
fact , all the three calculations give approximately the sarce 
u (x), except for a small region neaj: the inlet. 

Pressure gradient is shown in Figure 12. The three 
calculations are markedly different neam the inlet, hut 
slightly away from it the three curves merge into a single 
curve and they approach- a common asymptotic value. |dp/dx! 
is expected to he large near the inlet. With increasing 
number of terms, values of jdp/dxj predicted are indeed 
larger near the inlet, an-d show a sharper trend. The area 
under the dp/dx curve gives the pressure aaid is shown in 
Figure 12. 

The results of Morihara (1972) for R = 1000 are also 
shown in Figure 12. The streara\s/ise pressure gradient and 
as vjell the pressure are functions of both x and y coordi- 
nates in the above calculation. However, they are almost 
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independent of the transverse coord inc?te for large x. The 
points in Eigure 12 are those corresponding to = 0. The 
present results compare favourably vith the pressure gra- 
dient obtained by Morihara. Rapid decay of |dp/dx; near the 
inlet is approximately represented by the present calcula- 
tion. The pressure distribution of Morihara. is in good 
agreement ¥Jith the 5-terni calculation. 

The asymptote to the curve p(x) is given by equa,tion 
(3.53) j where is the pressure recovery coeffioient and 
is the value of (p + 3x) as x-^-oo . p^^ represents the 

value of ( p + 3x) at the point upto which the calculation is 
carried (see Equation (3.55) ). P^ is obtanned using 
Equaition (3.54). Table 33 - contanns the values of P^^ and 
P^ . The difference between the values of any two caJcula- 
tions is the area between two respective dp/dx curves. With 
the increasing number of terms j the prediction of dp/dx is 
expected to converge near the inlet with an associated 
reduction of the differences in the numerical values of 
the pressure recover 3 ?- coefficient. 


¥. 

^r1 

-P 

r 

2 

0.2187 

0.2437 

w 

0.2645 

0.2913 

5 

0.3245 

0.3505 


TIBIE 3 a 


PEESSUEE EECOYSEY COSillCISlT 
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In Table 31 the pressure recovery coefficient is 
compared \^ith tha,t given by other investigators. Collins &. 
Schowa.lter (1962 , 1963) vjho had earlier obtained a value 
0.338 have later reported 0.37. The calculation of 
Morihara (1972) 8.t E = 1000, gives a value of 0.365. In 
the first case of Wang & longwell (1964) there is a 


Eeference 

-P 

r 

Schlichting (1934) 

0.301 

Bodoia & Osterle (1961) 

0.338 

Collins &. Schov/alter (1962) 

o 

CO 

" (1963) 

0.37 

Vang & Longvvell (1964)"^ 

0.370^ 


0.352" 

Morihara (1972) 

0.365"^'’ 

Pr e se nt 5 -t e r m 

0.351 


TiEIE 3b COMI'lEISOE OE EEESSUEB IlECOVEEY COEIEICIMTS 


1. Values refer to case II, E=75> and are averages. See P.55 
2 . Inc lud e s up s t r e ara c o nt r ib ut ion. 

3. Excluding upstream contribution. 

4. Ecr E = 1000. 
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signif icejit vajriation of pressure across a iea:-dovmstream 
sectiorx, vvhich is attributed by the authors to the dis- 
continuity at the inlet. Also, an error has been committed 
in the a.pplication of the boundary condition an vorticity 
is assumed to be zero at the inlet. In the second case, 
pressure is more or less constant across any section suffi- 
ciently avjay from the inlet. Three values of pressure 
given at y = 0.1, 0.5 and 0.9 are averaged. Also, these 
values are relative to pressure at x = -00. ■ Therefore 
includes a contribution from the upstream of the inlet. 
Thus, the value for excess pressure drop obtained ’flith five 
terms is somewhere in the middle of the range of the values 
given by other ^'jorkers. ¥e , therefore , find that given 
by the present calculation agrees quite well with thnt 
obtained from calculations based on NSE for large R. 

The wall-vorticity o; = - u at the upper wall is 

+ y 

shown in Pigure 13. It is also proportional to skin 
friction. He suits due to liorihara for E =- 1000 are also 
shown in the figur-e. The wall-vorticit 3 ' is infinite at the 
inlet due to the form of the assumed initial velocity 
profile. The wall- vorticity compcires well with that 
obtained by Morihama except for a small distarxco neax- the 
inlet . 

Velocity rwrofiles at various x based on the 5-term 
calculation are shown in Figure 14- The 3 ^ show a small 
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waviness near the irxlet due to the representation of the 

given ve‘locity profile. 'Ihe vjaviness rapidly smoothens out 

and the velocity profile appears like a groviing -wall 

hound ary- layer together with a uniform central core, a.s 

assumed hy Schlichting (1934), (figure 14, see middle row). 

This comparison can he further checked hy examining the 

2 

growth of (u - 1) for this intermediate range (.004<x/ 
0.02), (Pigur e 15). The clear linear trend (cf . equation 
(9.69) of Schlichting, 1968) shows that the method brings 
out the qualitative features of che flow. The constant of 
proportionality obtained from the graph is 1.61 in compari- 
son with the Blasius constant 1.72. The bottom row of 
Figure 14 shows how the velocity prof iles evolve into 
Poise uille flow at distances away from the inlet. 

Pevelopiment length as defined in section 3.7.1 
is given in Table 4 atlong with those ohts.ined hy other 
?;orkers. The present values appear to he slightly less 
than those reported hy others except that of Schlichting. 
Development length based on the 5-term ca.lculat ion agrees 
with the value of Morihara (1972) within Ifo. 

Figure 16 comiTjares the 5-term values of u^ with 
those of Bodoia & Osterle (1961), and of YTeng and longwell 
(1964) for R = 75, and those of Morihara (1972) for R = 1000. 
(Some of the points of Wsng £ Dongwell shown in Figure 16 
are taken from Figure 3 of Vai Dyke, 1970). The results 
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Reference 

^99 

Present 2-term 

0.1642 

5-term 

0.1662 

5-term 

0.1679 

Schlichting (1934) 

0.160 

Hwang & Pan (1964) 

0.1688 

Morihara (1972) 

0.1692 

Bodoia & Osterele (1961) 

0.176 

Gillis &. Brandt (1964) 

0.1768 

Roidt &, Cess (1962) 

0.1816 

TABID 4 BEVELOPMT 

IDi'JGTH 


of Wang and Longvjell and of Morihara agree well with that of 
Bodoia and Osterle near the inlet. Hovjever , away from the 
inlet they agree more closely with the present calculation. 
Is mentioned earlier, the difference near the inlet could 
he due to the very smcill number of terms. Possibly 
agreement there could be improved by including more terms. 
Incidentally the limit process used to obtain the governing 
equation corresponds to the downstream limit. Therefore 
the present calculation could be at slight variance, near the 
inlet, with the results obts.ined for some large value of R. 

A similar agrument also holds for the calculation of 
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Bodoia & Osterle . But however one finds that the agreement 
of their calculation with that oased on BSE for large E is 
"better near the inlet than away from it. This difference 
could "be due to the accumulation of error as pointed out 
earlier . 

It should be iioted that the criterion of determining 
the initial condition leads to a good agreement slightly 
away from the inlet, although the predictions are not as 
accurate near the inlet. That is to say in some integra^l 
sense the initial condition in the present calculetion 
represents the flow well enough for even a few terms to give 
quite accurate predictions slightly away from the inlet. 
Inclusion of more terms may hamdlj' affect the calculation 
away from the inlet (figure 11), however, we expect 
imxDrovement in predictions near the inlet. 

'As the basic equalions in the present work are the 
same a^s those of Bodoia & Osterle (1961), coranarison with 
their work can be regarded as - check on the accuracy of the 
present method of integration. Also it might be argued that 
a systematic difference in rwo curves shown in figure 16 
might be due to e, shift in the initial condition. The 
forw/ard-difference formulation of Bodoia & Osterle could 
lead to significant accumulation of error. Tfe therefore 
regard the solution of the complete equations a:,t large 
Reynolds number to be a reliable guide in establishing our 
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method . Then the large x heheviour can he said to he well 
predicted. Since the prediction near the inlet can he 
improved using more terms in the calculations a good 
agreement for all x with calculations based on MSB at 
large E is possible in principle. 

It has been demonstrated in this chapter that the 
method predicts the entry flow in a straaght channel quite 
wells a.lthough the maximum number of terms used in the 
calculation is limited to five. IJot only the quantitative 
comparison of the centreline velocity, the pressure, the 
waBl-vort icity etc. are favour aiole, but also the calculated 
flow shows qualitative features like a. wall boundary-layer 
and a uniform, central core over a certa.in intermed ia^te 
range of x. Of course, the real test comes in the next 
chapter ?;hen we applj" the method to flows having 


r e c ir c ulat ion. 



CtlAPISR 5 

APPLICATIOK 20 SSPAPAIED PM PROBLEMS 


5.1 BTTROLIJCTIOR 

Ha'ving tested the method of calculation for the entrj^ 
flow prohlem in the pre^zious chapter, ?je give results in this 
chapter for separ'ated flow problems where many earlier method 
of calculation vmuld not be normally applicable. 

P3.0W problems in three types of geometrical configu- 
rations are considered (Pigure 1). The first type (figure 
la) is the flow in a chaiinel having symmetric sudden 
expansion. The second type of problem (figure 1b) is flow 
in an as3m3metric sudden expansion, vMiich can also be 
regarded as a baclaivard -facing step with a. top wall, or a 
base with a long splitter plate in a, chaimel. The third 
type of problem (figure 1c) deals with a base placed 
symmetricall^v between two parallel flat walls. 

We further consider two types of entij- condition in 
each of these problems. A parabolic velocitjr profile 
corresponds to a full 3 r developed entr 3 '" ccnd ition. A short 
inlet lea.ds to a miiform velocity profile. It is ^cossible 
to consider intermediate cases having anon-zero bound ar^^- 
layer thickness a little upstream of, say, the base. It 
is also possible to consider other geometrically different 
problems. However, we restrict ourselves to what seems to 
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"be a representative sample of the ahove three types of 
problems with one geonaetrical parameter h and with two 
simple hut extreme entry conditions. 

The number of terms N in (3.33) has to be determined 
for each calculation. ¥e can increase 11 till any further 
increase does not affect the results significantly. Alter- 
natively, we can fix K and find the range of the geometri- 
cal parameter in T^hich results are acceptable. 

Consider now the flow in a symmetric sudden expan- 
sion with h = 0.5 when the entry profile is parabolic. ITS 
solution for values of R upto 160 is given by Macagno and 
Hung (1966). Their figure 7 shows that an asjmiptotic trend 
has been attained atid therefore their results can be used 
for comiparison with the present calculation. Morihara. (1972) 
has also obtained ITS solution for this problem for values 
of R upto 20 (see figure 5). The rjresent calculation is 
started with a few terms and later more terms axe included. 
We shall see that five terms are sufticient to indicate 
convergence . Therefore , the number of terms in subsequent 
calculations IS restricted to five. 

At this stage it may be appropriate to briefly 
summarise the important steps in the development of the 
present work. A limit equation is obtained for large R for 
a class of problems which may have recirculating regions. 
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A method of flov; ca.lculation is then dove lotted vfnich leads 
to a set of ordinary differential equations, lije upstream 
condition is used directly to obtain the required initial 
condition. The concern of the present chapter is to test 
the consequences of the basic inputs, that is the calcula- 
tion rrethod in its final form when the flov^s have recircu- 
lating regions. It is applied to a variety of situations 
and comparison with available results is made wherever 
possible. Also some of the results are believed to be new. 

5.2 CHAlWiEL A SYIaI\3ETRIC SbEDMI EXPAl^SION 

Consider the flow in a channel harming a symmetric 
sudden expansion (Figure 1a). The expansion ratio, which 
is defined as the ratio of the dovmstream charnel width to 
that upstream, is given by 1/h. 

5.2.1 HOT I AL COMFIT lOli 


in the case of the parabolic velocity profile at the 
entry, the initial flow is given by: 


u^ = 3/(2h)(1 - iy/hf) 


= 0 


j 0 s 

, h p |y ; y 1 . 


(5.1) 


The initial condition is then obtained to be 


m 


a ^ = -5 , h b - t 

mO m ffi 


/ (s F) 


y-=h 


m 


(5.2 
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where 

1 

• ( 5 . 3 ) 

0 .. 

When the inlet velocity profile is -uniform, that is 
= 1/h , 0<|y[<h , 

= 0 , h^;y!^1 . (5.4) 

the initial condition for a is fo-und to he 

m 

^mo = - C 5) • (5.5) 

We first discuss the case for which NS sol-ation has 
been obtained by Hung and Macagno (1966), Other cases are 
discussed later. 

5,2.2 ESSUIT'S FOE h = 0.5 vfITH PARABOIIC ENxEI PROFIIE . 

The assumed parabolic entry profile given by (5.1), 
is shown in Figure 17 with its approximation using 2, 3 
and 5 terms. Clearly, the given velocit;/ profile is better 
8.pproximated with increasing number of terms. 3- end 5-term 
representations are closer to each other. 

Streamlines based on the 3-term calculation are 
shown in Figure 18. Also streamlines obtained by Hung and 
Macagno (1966) for E = 46.6 are replotted with the present 
variables in the same f ig-ure . It is seen that the stream- 
line pattern obta.ined by the present method is in excellent 
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qualitative agreement with that obtained by Himg and 
Macagno . The numerical details of the flow are compared in 
Table 5. Values given in the first row are obtained from 
the graphs given by Hung and Macagno. 



Re attachment 

Centre of 

Recirculgt ine 


point 

eddy 

me.ss flow 



z 



r 

0 

■ e~1 

Hung and 

Macagno 0.066 

0 . 01 3 

0.052 

Pre sent 

3-term 0.064 

0.014 

0 . 045 

TAEIE 

5 SYMMETRIC SUDPEH EXPANSION , h = 

= 0.5 

JL ij. X S 

favourable comparison 

of the preser 

it ca,lcula- 


tion with tha.t of Hung and Macagno (1966) which is based on 
NSE, supports the following views. Equation (2.5) is 
appropriate an the limit equaition for large E for the 'far' 
region, and the streanr^/ise length scale is of order R. The 
method devised to -solve equation (2.5) is applicable even 
in the presence of recirculating flow. Also the hypothesis 
used to obtain the initial condition to integrate the 
equation aippeans to give acceptable results. 

The centrelirig velocity u is shown in Figure 19. 

c 

3-term and 5-term calculations give almost the same u 

c 

except very close to the expansion. Although a very small 
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n'um'ber of forms are used, the results c?n oe srid tc be 
indicative of a convergent trend. Velocity profiles a.t 
various streaaw.!ise locations based on the 3-tGrm calcula- 
tion are shown in figure 20. The velocity profile e.t 
X = 0.06 is apijroximately that at the rea-ttachment point. 

figure 21 shows the pressure distribution obtained 
using 2, 3 and 5 terms. The point of re attachment is 

marked on each of the curves. At first the pressure increa 
ses followed hy a continual decrease a,s required by the 
Poiseuille flow. It is seen that the mancimurn pressure 
occurs after the reattachment . The numerical values of 
^99 given laiter in Table Sc vjith those for 

other values of h. 

5.2.3 EBSUITS POPl OTHER VALUES OP h 

Numerical calculation with two types of entry 
profiles, the pera.bolic ai'id the uniform, are carried out 
for various values of h given in Tahle 6. a zero indi- 
cates the occurrence of the singular it 3 ’' and the terminstio 
of the calculation (Sections 3.5 end 3.6). h = 1 corres- 
ponds to the entry flow problem. 2-term calculations 
ale also carried out for some intermediate values of h 
in addition to those given in Table 6. Table 7 gives 
some typical initial vsl-ues of 




N h = 

2 

3 

5 

2 1 

3 1 

5 1 


0.95 0.9 0.75 0.5 0.25 0.1 

Parabolic entry profile 

1 11 11 1 

1 11 10 0 

1 11 10 0 

Uniform entry profile 

1 11 11 1 

.1 1 1 1 1 0 
1 11 10 0 


■TABLE 6 CALCUIATIOMS PEFiFOPlIEL POR 
SMBTRIC SULLEU SXPMSIOK 


i s D 5^ 0 

G t D = 0 for some x (see Section 3.5). 
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^•10 ®20 ^30 ^40 ®-50 
Parabolic entry profile 


0.9 

.1322 

-.0160 

.0330 

-.0073 

m 

O 

0.5 

1 .1668 

-.0200 

-.2516 

.0228 

.0650 

0.1 

3.8088 

1 .0040 

3.1894 

.9347 

2 .3766 



Uniform entrj' ; 

profile 


0.9 

-.3751 

.0443 

-. 0865 

.0176 

.0348 

0.5 

.6744 

-.1758 

.0059 

.1124 

-.1807 

0.1 

3.7262 

.9442 

2.7992 

.7358 

1 .5681 


7 imiAL COE'PIilOi; POR 
SYISIEIRIC SURRM EZPMgioy 
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Foi; small values of h the initial flow is like a 
strong jet, and therefore the higher modes are expected to 
ho significant . On the other hand , for h close to one 
the deteiils of the flovi/ will not he hrought out accurately 
unless higher modes are included in the calculation. For 
b = 0.1, the order of. magnitude of (see Table 7 ) 
suggests that more terms should he considered in the expan- 
sion. On the other hand for h = 0.9, the order of magnitude 
of suggests that inclusion of more terms may not he 
necesseny. liowever , it turns out that the small 'hack-step' 
requires higher modes to he included for the calculation 
of flow details like recirculation. In the case of h = 0.5, 
which is in the middle of the range of h, the inclusion 
of the fourth and the fifth modes affect the calculation 
only slightly. 


For h = 0.5 the assumed uniform entry Vc;locit 3 ?- 
profile ( 5 . 4 ) with its a;,pproxima.tions is shown in Figure 22. 


The assumed velocity profile has 


disc ont inuit y 


5. 


The approximation is considerahly improved when higher modes 
are included, iniother exeinple where the assumed velocity 
profile is discontinuous is the entry flow in a channel 
considered in the previous chapter. However, there the 


d is c ont inuit y is 


af the wall, ’whereas in the present case it 


is located within trie interval (0,1) 
approx ima,t ion is improved ‘with the L 
(see also Figure 10). This suggests 


j.n Doth cases the 
elusion of higher mode 


that discontinuous 


o 



69 


velocity profiles can "be expanded into eigenfunctions of 

(3 -6). figure 22 also shovus n obtained vith 2, 3 and 5 

c 

terms. The three calculations show a consistent trend. 

Pigure 23 shows u for three values of h based on 

c 

5”term calculations with both types of entry conditions. 

Por the paf abolic entry profile the- figure shows a monoto- 
nic decrease of to the Poiseuille flow- value. An 
interesting behaviour is observed for h = 0.75 with uniform 
entry profile. Pirst decreases and then increa,ses to 
attain the Poiseuille flow value. The incoming flow, due 
to its mixing with the ’dead' fluid, first decelerates, 
later as the flow develops, u increases as u (0) u (c.-). 
The above qualitative trend is not shown by the 2-tera 
calculation. We therefore see that the inclusing of higher 
modes in the calculation helps in realising the flov? picture 
more accurately and thus a, better quantitative prediction 
can be obtanned. Por h = 0.9 the 'dip' in u^^ is however 
absent. Possibly the 'dip' is small and -may require more 
terms in the calculation. Pigure 24 shows the pressure 
distribution based on 5-term c ale ula-t ions for the above 
values of h. While each of the curve for the parabolic 
entry profile show a positive pressure recovery coefficient 
positive as well as negative occur for uniform entry 
profile (see Tables 8a, b) . Points of reattachufint | x^ have 
been marked in the figure. Pollowing x^. a maximum. 
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h 

N 


^99 

X 

r 

0.95 

2 

.0445 

.0712 

- 


3 

.0521 

.0784 

- 


5 

.0585 

.0706 

- 

0.9 

2 

,1033 

.1157 

- 


3 

.1199 

. 1 1 48 

- 


5 

.1360 

.1140 

- 

0.75 

2 

.3892 

.1876 

. 0030 


3 

.4292 

.1862 

.0080 


5 

.4622 

.1860 

O 

GY 

O 

O 

0.5 

2 

1 .320 

.2839 

.0659 


3 

1.234 

.2779 

• 

O 

00 


5 

1.390 

.2775 

.0656 

0.25 

2 

2.022 

.3594 

.1431 

0.1 

2 

1 .838 

.3715 

.1558 

TABIE- 

8a SmmjiRY OB 
Pi^lBOLIC ENTRY 

EESUITS WITE 
PROEIIE 




h 

K 

P 

r 

^99 

X 

r 

• 1 

2 

-.2437 

.1642 


(Entry 

flow) 

3 

-.2913 

. 1 662 

- 


5 

-.3505 

.1679 

- 

0.95 

2 

-.2604 

.1673 



3 

-.3147 

.1693 

- 


c; 

-.3699 

.1711 

- 

0.9 

2 

-.2700 

.1696 

- 


3 

-.3179 

.1714 

- 


5 

-.3570 

.1726 

- 

0.75 

2 

-.1590 

.1559 

- 


3 

-.1096 

.1571 

- 


5 

-.0955 

.1673 

• 0007 

0.5 

2 

.6397 

.1495 

• 0219 


3 

.6889 

.1438 

.0234 


5 

.4970 

.1386 

.0220 

0.25 

2 

2.028 

.3432 

.1265 


3 

2 .285 

.3406 

.1253 

0.1 

2 

1 .884 

.3708 

.1551 

TABEE 

81 

SLloMi-RI OP 

EESUITS 

WITH 


Uiai'OPul EKTEY EROPIIE 
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occurs, as though there is a region of approximately uni- 
form pressure after the re attachment . A similar trend is 
observed in other calculations to be described later. 

Tables 8a and 8b give the va.lu 0 s of p , and x for 

pa,rabolic and uniform entry profiles respectively. 

5.2.4 PEESSUTiE REiCOVSEI COEEPICIEiiT 

It would be of interest to sec how the calculated 
value of P^ compares with an approximate integra.l ceJcula- 
tion. One such approximation is given here. Plow at a 
downstrea.m section is given bj' Poiseuille flow, and at x = 0 
it isal^own flow. The constant pressure gradient of -3 of 
the fully developed flow and any direct viscous effect are 
disregarded. Pressure P^ acting on the downstream section 
is therefore due to the change in moiTientijm fliuc. 

Y^e now calculate P with the above approximation. 

r 

When the entering flow is parabolic ( 5 . 1)5 momentum balance 
between the sections at the expansion and at some distance 
downstream gives 

= 6/5 (l/h - 1) . (5.6) 

Y<hen the entering flow is uniform (5-3) s ’"Q gel' 

P = 1/h - 6/5 . (5.7) 

X 

Approximate calculations have been earlier based on 
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■uniform -velocity profiles which give 

= 1/ii - 1 . (5.8) 

It should he noted that in (5.7) could he positive or 
negative depending on the value of h, whereas (5.8) always 
predicts a positive value for P^. 

The approximations on v;hich the above formulae are 
based become increasingly inaccurate as h becomes very smal 
or the expansion ratio becomes ver5'- lamge . They can be 
expected to provide an estimate for moderefe expansions. 

The calculated values of P are compared with the 

r 

approximate values given by equations (5.6,7) in Pigure 25 
To avoid overlapping some of the points in Pig-ur-e 25 (end 
also in Figure 26) have been staggered. Por a small 
'back-step' i.e. h close to one, with the uniform entry 
profile, P^ is approximately equal to that for the entry 
flow problem (h .= 1). Por the parabolic entry veloc it y 
profile P^ goes to zero as h goes to one . It is seen that 
for the expansion rafio upto about two , equations (5.8, 7), 
sho-wn by solid lines, predict P^ rather ■well. As h goes 
to zero these equations are singular and therefore cai'inot 
be expected to be valid for large expansion ratios. 2-term 
calculation predicts an. approximately uniform level of P^ 
for large expansion ratios. P^ is expected to be reduced 
considerably as the losses due to the expansion will be 
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significant for large expansion ratios. It is sioggested 
that for the- range of the expatision rajtio greater than, say, 
three , may he approximately taken to he a constant and 
equal to two- irrespective of the entry condition. 

5.2.5 KEAflACmaEM' Ai^D DEYEIOBVrEl'jI lEhC-THS 

The reattachment points x and the develonment 

r 

length Xg^ are shown in Pigure 26. Por h = 0.9 ana 0.95 the 
reattachment is not obtained with upto five terms. As 
noted earlier, more terms will he needed for smaller ’hack- 
step' , to predict accuratelj/ the details such as reattach- 
ment. Por the ratio of expansion equal to three, Durst 
et . al- (1974) have calcula.ted the flow using NSE. Although 
the upstream condition used in their calculation is derived 
from the experiment, it is believed that the entry velocity 
profile is approximatly parabolic. I‘he point of reattach- 
ment as read from their Pigure 6 is shown. Although R in 
their calculation is only moderately large, x^ obtained by 
them compares well with the present result. Generally x^ 
with the uniform entry profile is smaller compared to that 
with the parabolic entry profile. It is to be expected that 
the flo¥; reattaches earlier in the former case as the 
mixing layer is stronger thaXi in the latter case, x^p 
obviously approaches zero with parabolic entry profile as h 
goes to one. However with the uniform entri' profile the above 
limit corresponds to the entry flow problem. 
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^.3 CliANl^I m ASmffiTEIC SUDDEN S}:riNSION 

The flow in a one-sided expansion of a channel shown 
in figure 11 is novj considered . The primary aim of this 
application is to demonstrate the use of antisymmetric 
e igenfunct ions . 

5.3.1 ffllTIAl CONDITION 


The downstream channel ^i/idth is two. The step 

height is denoted hy h. Ylhen the entry velocitj^ profile- is 

parabolic the initial condition for a becomes 

m 


m 


®'m0 


1.5 X (-1) 
(1 - h/2)^G 


m 


(1 - h/2) - 0^ I 

m EiJ y = 1 - h 5 


(5.9) 


0 £ h ^ 1 


where G is given by (5.3). In the case of a uniform entry 


m 


profile, the initial condition for a is obtained to be 


^(- 1 ) + (- 1 ) 


^•mO “ 


m 


’ m 


(1 - h) ) 


(2 - h) G, 


, 0 -4h £ 1 . (5.10) 


m 


Equations (5.9) and (5.10) can be modified when 1^ h £ 2 


5.3.2 EESUITS 

Table 9 gives various values of h for which the 
calculations are carried out. With pamabolic entry profile, 
the given flow for small h differs only slightly from the 
fully developed flow (see Table 10). The small perturba- 
tion results given in Section 3.3.4 can be used, for h = 0, 
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u 

h = 0.05 

0.1 0.25 0.5 

0.75 

1 .0 



Parabolic 

entry profile 



2 

1 

1 

1 1 

0 

1 

3 

1 

1 

1 1 

1 

1 

5 

1 

1 

1 1 

1 

0 



Uniform 

entry profile 



2 

1 

1 

1 1' 

1 

1 

3 

1 

1 

1 1 

1 

0 

5 

1 

1 

1 1 

1 

0 


TABIE 9 CALCULATIONS PEEPOPJ/lE-D 
ASm®TRIC SUDDEN EXPAI'jSION 

FOR 



1,0; refer 

to Table 

6. 
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h 

^-10 

^20 

a^O 

o 

\ 

^50 



Parabolic 

entry profile 



1 .0 

-.8595 

-1 .0544 

-.0256 

.0727 

-.0595 

0.5 

-.3877 

.3129 

-.1066 

.0080 

.0035 

Q.1 

-.0108 

.0591 

-.0067 

.0071 

-.0069 



Uniform e 

ntry profile 



1 .0 

-.6062 

-.3340 

.3203 

. 0401 

-.3017 

0.5 

.1670 

.0021 

-.1056 

.0319 

-.1524 

0.1 

.0659 

-.3535 

.0404 

.0421 

.0405 


TABU3 10 HOTIAI COUBITIONS BOR 
ASYIvffiffiSRIC SIIDDM EXPANSION 
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the problem is trivial. In the case of uniform entry 
profile, h = 0 corresponds to the entry flow problem. 

Initial conditions obtained for some typical values of h 
are given in Table 10. ' 

The given para.bolic entry profile with its approxi- 
mations is shown in ligune 27 for two values of h. \;ith 
increasing number of terms the approximation is seen to 
improve. In the figure instead of showing each of the 
three approximations (2, 3 and 5 terms), only two have beer- 
shown. The omitted, one lies close to the next higher 
approximation. 

When the entry profile is uniform, its approximations 
are shown in Figure 28. Sven 5 terms provide a poor appro- 
ximation in this case. In fact for some values of h it 
appears that the total flow is less in the upper half of the 
channel than in the lower half, whereas the opposite is 
true for the given profile. However, a convergence of the 
oscillatory type cannot be ruled out. Here, the numerical 
results for only the parabolic entry profile are discussed. 

The asymmetry of the flow is shown in a simple way 
by considering the skin friction at the two walls. The 
values at the upper and the lower walls are denoted by 
suffix + and - respectively. The skin friction coefficient 
is also related to the wall-vorticity 60 by 



79 




R = + 0^ 


+ 


+ 


i c^R 3 .t the tvv’o ¥;alls aie shown in Riginre 29 for h=1 . As 
the flow develops the 8.synimetry of the flow reduces and vaiii- 
sb.es ult imate ly , Compared to the 2— term calculation, the 
3-term calculation shows more asymmetry/ a small recircula“ 
ting region is obtained on the lower wall with three terms. 


The point of reattachment is taken to be the point ■where the 
skinfriction curve meets the x axis going upwamd (see also 
Figure 30). 


Velocity profiles at various strearawise locations 
based on the 3-term calculation are sho'^vn in Figure 30. The 
gradual evolution of the flow into a Poiseuille flov^ is seen 
from the f ig'ure . Pressure distribution is shown in Figure 3'I. 
Predictions by the 2~term and the 3-term calculations are 
qualitatively similar. At the begirning the pressure drop is 
large followed by a region over which the pressure is appro- 
ximately uniform., and then the pressure decreases further to 
its asymptotic trend. It is interesting to note that the re- 
attachment as predicted by the 3-term calculation marked by x^ 
in the figure , is followed by a region of approximately uni- 
form pressure. Pressure distribution is shown in Fig-ore 32 
for sotiie typical values of h. Numerical results are summari- 
sed in Table 11 ?;hicli gives P^ , x^^ and x^ for various valuJ" 
of h. 

It should be noted that the eigenvalues are numbered 
in the increasing order, and. the first and the third 
eigenfunctions are antisymmetric. 
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h 

N 


I 

I 

I a \ 

X 

1: 

0.05 

2 

.0147 

m 

o 

00 

- 


3 

.0147 

.0318 

- 


5 

.0209 

.0318 

... 

0.1 

2 

.0333 

.0751 

- 


3 

T“ 

CM 

O 

• 

.0732 



5 

.0455 

.0726 

- 

0.25 

2 

.1078 

.1428 

- 


3 

.1017 

.1361 

- 


5 

.1388 

.1319 

- 

0.5 

2 

.3114 

.2101 

- 


3 

.2553 

.1892 

.0094 


r~ 

0 

.3437 

.1883 

.0145 

0.75 

3 

.4104 

.2463 

.0409 


5 

.4252 

.2365 

.0421 

1 .0 

2 

.1617 

.2857 

- 


3 

.1016 

.2799 

.0617 


TABLE 


11 ASYIMBTRIC SUDDEN EXPANSION 
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5.4 BASE M A CHASEEL 


Consider a "base ha'ving thickness 2h, placed symmetri- 
cally in a channel as shown in Eigure 1c. Ihe flow he hind 
a thick long flat plate with a bluff end placed in a v\?ind— 
tunnel can he idealized to lead to the present problem, 
although it is recognised that only a few features of the 
base flow are modelled in the present problem. 

5.4.1 miflAI CONE n ION 


Two types of entry conditions are chosen as before. ' 
Eor the parabolic velocity profile at the entry, the initial 
condition for is obtained to be 


6 


mO 


(1 - 


+ 2/(1 - h) 


. (5.11) 


-'y 


h 


In the case of the ruiiform entry profile the initial condi- 
tion for a turns out to be 
m 


1 


a 


mO 


(1 - h)G 
^ m 


f)"(h) - e"(1) 


( 5 . 12 ) 


5.4.2 RESULTS 

Table 12 shows the range of the calculations. Table 
13 gives initial conditions for some typical cases. It 
should be noted here that as h 0 the base shrinks to a 
thin plate . Eor the parabolic entry profile the problem 



W h = .05 .1 .25 .5 .75 

PaialDolic entry profile 
2 11111 

3 1110 1 

5 1110 0 


Uniform entry profile 

2 
3 
5 

TABLE 12 CALC UIAT 1017 S PERFOEIED 
FOR THE BASE 117 A CHANIEL 


1 11 1 1 
1 1111 
1 0 1 1 d) 


1 , 0 


; Refer to TaLle 6. 
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h 

^■10 

^20 ^30 

^40 

^"50 



Parabolic entry profiP 



.1 

-1 .4384 

-.0690 -.0769 

.0016 

.0250 

.25 

-1 .7555 

.0738 .2185 

.0296 

-.0377 

.5 

-1 .3457 

.3270 -.5916 

.0205 

.1490 



Uniform entry profile 



.1 

-.7851 

.0604 -.4038 

.0603 

-.2316 

.25 

-1.3608 

-.0630 -.0646 

.1101 

.1305 

.5 

-1 .3424 

.2559 -.1729 

.0737 

.0774 


TABIE 

13 fflITIAL COEUITIOU POP 
BASE B! A GHANBEI 

THE: 
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considered is then like a flow behind a trailing edge of a 
flat plate hawing finite, equal and opposite vorticity on 
the tv\/o sides, vi^hen the entering flow is uniform, there 
are two vortex sheers of equal and opposite strength. Ihe 
e igenf line t ions which are used in the approximation have zero 
vorticity at the centreline. Therefore for h = 0, the 
present approximation is non-uniform in a, neighhoiirhood of 
y = 0. We observe that as h-> 0, expression (5.12) 
reduces to (4.2) which corresponds to the entry flow problem. 
However for h 7^ 0, the approximation is expected to have 
the desired features of the given entry profile such as 
u = 0 at y =: 0. 

Figure 33 shows the entry condition as approximated 
by 2 5 3 and 5 terms for h = 0.25 with the parabolic entry 
profile. Table 14 compares the centreline velocity at the 
entry u and the associated momentum flux M given by the 

three approximations with the actual value. Once again we 
see that the approximation is improved with the inclusion oz 


2-term 

3 -terra 

5 —“t 0 T Til 

Actual 

0 

0 

[ 

« 

00 

.0368 

.0287 

0 

1.6139 

1 .5978 

1 .5984 

1 .6 

TABIE 14 REPEESM' 
EOE THE 

TITIOH OE 
BASE IE A 

THE lE'ITIAL 
CHilEEl 

mourn 
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higher modes. 

The approx imai ions io uniforni eniry condition for 
h =0.1 and 0.5 are shown in figure 34. Although any marked 
improvement in the approximation with the inclusion of 
higher modes is not evident for h = 0.5 , it is cleamly seen 
for h = 0.1. 

The streamline pattern "based on the 5-term calcu- 
lation is shown in Figure 35 for h = 0.25 with parabolic 
entry profile . A closed "buh'ble "behind the "base is o"btained . 
Near the "base the streamline pattern shows slight waviness 
which is possi"bly due to the representation of the initial 
condition. This waviness decays rapidly. The development 
of streamwise velocity profiles "based on the 5-term calcu- 
lation is shown in Figure 36. The centreline velocity is 
shown in Figure 37. The initially decreasing u^ due to 
entrainment is absent in the 2-term calculation , whereas 
3- and 5-term calculations give a better picture of the 
recirculating flow behind the base. Entrainment of the , v 
'dead' fluid follovved by detrainment is seen from the 
figure. Although the size of the recirculating region is 
approximately the same for 3- and 5-term calculations, the 
eddy is stronger and its centre lies closer to the base in 
the 5-term calculation. Table 15 compares the three calcula- 
tions. As the 2-term calculation does not lead to closed 
streamlines, some of the entries in Table 15 are indicated by- . 
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2-term 

3-terffi 

5-term 

Point of reattachment, x 

’ r 

.0086 

.0086 

.007 

Maximum half -width of the 
recirculating region 

.263 

.205 

.264 

location of the centre 
of the eddy, 




X 

^ e 

- 

.0042 

.0024 

^e 

- 

.108 

.146 

Total recirculation, - *4^ 

-- 

.007 

.021 

Maximum reversed flow 
velocit 3 ?- at the centreline 

-.1817 

-.0918 

-.2198 

Development length, x^q 

.213 

.206 

.201 

Pressure recovery 

-.492 

-.214 

-.116 


coefficient, p 

’ r 


TABLE 15 BASE III A CHABBEL (h = 0.25) 
?/rT'H PAEABOLIC ENTRY PROEIIE 
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Th6 prBSSuic© .d isti r ilDiiij ion is shown in Figure 38. The poini 
of re attachment is also marked on the figure . The 
figure suggests that there is a tendency to form a region 
of approximately uniform pressure following the point of 
r e at t ac hmen t . 

Figure 59 axid 40 show centreline velocity and pressu 
distributions respectively when the entering flow is uni- 
form. For h = 0,5 and 0.25, u^ decreases initially, 
however, no reversed flow is obtained. This evidently is 
due to the error in the approximation of the entry profile. 
It is interesting to note that for h = 0.5 the calculation 
shows a region of uniform pressure v^hich in the fewer terms 
calculation is absent. Tables 16a and 16b give , x^^ and 
x^ for parabolic and uniform entry profiles respectively. 
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ii 

N 

P 

r 

^99 

X 

r 

75 

2 

-.4835 

.1846 



3 

-.9255 

.2117 

- 

5 

2 

-.9233 

.2248 

.0216 

25 

2 

-.4921 

.2130 

• 0086 


3 

-.2138 

.2057 

.0086 


5 

-.1155 

.2007 

.0070 

1 

2 

-.1621 

.1955 

.0001 


J 

-.2262 

.1947 

.0005 


5 

-.1229 

.1948 

.0006 

05 

2 

-.1596 

.1944 

- 


3 

-.2561 

.1914 

- 


5 

-.1729 

.1916 

- 


TABES 16a 

BASE IN 

A CHANNEL WITH 



PAEilBOLIC ENTRY PROEIIE 
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h 

N 


^99 

X 

T 

.15 

2 

-.5190 

.1899 

- 


3 

-.9782 

.2157 

- 

.5 

2 

-.8727 

.2253 

- 


3 

-.9321 

.2264 

.0173 


5 

-.3219 

.2244 

- 

.25 

2 

-.1971 

.1970 

- 


3 

-.2360 

.I960 

- 


5 

-.2731 

.1901 

- 

.1 

2 

-.1789 

.1872 



3 

-.4337 

.1822 

- 

.05 

2 

-.2176 

.1795 

- 


3 

-.3854 

.1792 

- 


5 

-.4331 

.1765 

- 


TABLE 16L BASE IlM A CHAIBEL WITH 
UKflEOEIiI EKTRI PROFILE 



EPILOGUE 


EXIEEUAL SEPARAIEL I'LO'f 


Ihe internal separated flow provides a snitaLle 
grotind for the testing of limit equations for L^ ^ 1. le 
then find that for 1, two principal limits are possi- 

i 

hie and they correspond to 1,, -v' 1 and H. While we 8.re 

familiar with the first limit process vjhich leads to the 
Euler equation, the second limit process which is rather 
new is examined in this thesis. 

The limit equation (2.4) corresponding to the limit 
1 R is welllaiown and is associated with the classical 

X 

hound ar y-layer . We test the applicahilit}’- of this equation 
to a class of separated flows. Good comparison with solu- 
tions using complete USE supports its validity. 

External separated flows differ in one important 
respect from the internal separated flo’ws. The interaction 
between the outer irrotational flow and the separated 
region is an essential characteristic of the external 
separated flovi and poses considerable difficulty in the 
analysis. As it stands, it is difficult to ascertain 
apriori the transverse length scale of the separated wake. 
In fact some drastically different models of the limit flow 
have been proposed. However, if YJe assume that the tr,ans- 
verse length scale of the separated wake is of order one, 
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some of the conclusions of this v^orlc such as the appropri 
ateness of the limit equation (2.4) would be applicable. 

The experimental results gi^'/en b^r fcrivos and his 
colleagues (I964j 1965, 1968) for flow past a circular 
cylinder suggest 1 end E. The limit equation 

(2.4) therefore should be applicable. In this limit the 
■viscous anb the inertial forces are of seme order. This 
nature of the flow is different from the other shear 
driven recirculating regions where the limiting flow is 
expected to be governed by the Euler equation. 

The present work may give an insight into one 
element of the external separated flow at large Reynolds 
n'umber whose globa.l picture continues to fascinate many. 
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APFMDIX 1 


SOLUTION .OF THE EIGENVALUE EROELEM 


1 . The eigenvalue proLlem ohtained in Chapter 3 is 



A 1 3/2(1 - 



( 11 - 1 ) 



0 at y = + 1 



Calculation of odd and even eigenvalues proceed separately. 
The region of integration is restricted to Oi, y = 1 . The 
■boundary condition imposed at j = 0 is therefore as follovt?s 

Odd ; = n” =0 at y = 0 , (A.1.3) 

U 

Even : -p = eh = o at y = 0 . (H .4) 

For a given A tvo lineanly independent odd solutions 
of (a 1.1) can os obtained. One such solution is 



(A1 .5) 


Another solution is ootained on integrating (A1.1) with the 
initial conditions 



A' 1" ^ , 

at y = 0 



A general solution which is odd can then "be vaitten as 
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,(1 ; A) = 0 , 


(A1 .8) 


Thus the zeros of the function givef the odd 

eigenvalues. 

Similarly for determining even eigenvalues, 'we con- 
sider tvjo linearly independent even solutions of (11.1). A 
pair of such solutions is obtained on integrating (A1.1) 
v*ith the initial conditions 




-h 

■V3 

il 

0 

C V 
v_>i' ^ - 

= 

1 

— nci — r 

9 S- • ^ 

^ at 

y = 

0 

(A1 

,■ ", \ 

• t; 



4^.'. 

I 

' L 

= =:qT 

•4 t'4 

= 

0 

L nt 

at 

\r ZZ 

0 . 

(A1. 

10) 

If 

a non-tr 

ivial linear 

c 

o.mbination of 


and 

; 

s 

atisf 

le s 












(A1 

.2) 

, the 

n 











( ^-'3 

■ 4 3 

■ e 

7= 

= 0 . 

1 




(A1 . 

11 ) 

The 

le 

ft si 

de of (A1 

.11 

) 

is a function 

of 

X . 

Its 

zero , 


the 

ref 

ore , 

aives the 

even 

, eigenvalues. 








The 

values of 

4- In 
u ll 

e 

functions on t 

116 

left 

side 

of 


(A1 

.8) 

and 

(A1.11) can 

easily be obtsdned 

for 

a g ive 

n 

« 

In 

the 

integration of the 

equation (A1 . 

1) 

t he 

fourth 

ord e 

T 


Runge-Xutta-&ills method v’ith a step length 


of 0 . C5 


use a . 


Using the ilevrton-Raphson method zeros of the above fi'nction& 
are obtained . 


The eigenfunctions are scaled in the 


.CllOV'illl 



y = 0 
y = 0 


9 


96 


Odd i = 1 at 

, I 

Even s cip = 1 at 


(A1 .12) 
(A1 .15) 


Solution of the adjoint problem is obtained in a 
similar way. Note that 6= constant is a solution of the 
adjoint equation (3.23). 'The scaling used is similar to 
(A1.12, 13). 

Odd : 0 = 1 at y = 0 ^ (A1 .14) 

Even ; 0 = 1 at 3 ^ = 0 . (A1.1'~ 

Eng e invalues are given in Tables A1 end A2 . Each 
time a calculation is to be performed, it is not necessary 
to go through the process of recalculating the eigenvalues- 
They can. be given as data in the program and the e igenfunc't ■■ 
ions can be generated as solutions of an initial value 
problem. Three of the required four conditions at y = 0 
are given bj^ (A1.3, 12) or (a 1.4, 13). The fourth conditicr 
which can go as data in the program, is also given in 
Tables A1 and A2 . 


m 

1 14.4535 

2 48.8857 

3 104.572 




( 0 ) 


-5.82126 -8.45962 
5.70202 150.811 

-10.3516 -60.7067 


TABLE A1 


ODD EIGEIJVALLES 



lTi 

A 

IB 

c-D (0) 

' 171 

V 

’(0) 

1 

18.8146 

-18.0682 

-22. 

4208 

2 

57.5463 

-194.425 

-90. 

9318 


117.486 

-116.047 

-167. 

001 

L 

199.103 

17 G9 0 1 

^ 9> — • ■ 5 ^ 


759 

1— 

j 

303,073 

-312.985 

-441. 

431 


TABLE a2 EVSS EIGBITAIuES 
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11 =: 

1 

2 

3 

4 

5 

m = 1 

1 

.0003 

.0010 

.0014 

.0020 

2 

. 0009 

1 

.0058 

.0034 

.0129 

J 

,0009 

.0018 

1 

.0076 

.0086 

4 

.0030 

. 0058 

.0193 

1 

.0614 

5 

.0030 

.0036 

. 01 1 6 

.0192 

1 


TABIE A3 

A 0111:01 

OH BIOSlHOGOHALIiy 



adjoint proolem iDeconies rather large compared with distance 
o'^/er which significant ■variations occur when higher modes 
are considered. If we decrease the step size, we can expect 
this error to reduce signif icant ly. However we acceijt 
these calculated e igenfunccions and efforts to obtain them 
more accurarel;/ were not made. Inc id ent 3.11y this exercise 
confirms the development of the biorthogonal set. 


APIENDIX 2 


ON THS IHASE PLAICE 


The differentia,! equation (3.36) is written as 


(1 - ^l22®-2^' 




+ C2 ^ 2 ®' 2^®-1 ^ " ^ 221.®':1 ^ 222®2 


)a,-,- -\a^.(-n2.2) 
cl 


'i‘Ne above equations can be solved lor a^ and a.^ so give 


Cv^i — iJr^l Ij 5 


Vine re u? D. ana 


:e functions of and ^iven telow 


^^ 11^221 '^ 121 *^ 211^^1 


'^ 121 *^ 211^^1 ^ ^'^ 11l'^222 ■'■ '^ 112'^221 


*^ 122*^211 ^®' 1®'2 ^'^ 112*^222 " *^ 122 ^ 211^^2 


I 1 i 


222^^2 * 


(A2.4) 


0222 n'^'iii 


'122‘^2 


(•'^2^111 ^i*^212^^1®'2 


(12 J 


ex pr e s s i ons are q ud ar at ic , t hat i s t he y ar 9 o f t rp 
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(.42.7) 


The numexicsil Y&lues of the coefficients, /'s of D, and 


Dp ais given in Ts'cle A4 



jj 

^1 

JJo 

c. 

■ 1 

.1796 

-5 .546 

-1 .267 

/7 

■'2 

-.1793 

74.79 

-22.23 

'-'3 

-3.702 

-611.3 

-203.6 

.<-4 

.6947 

-18.81 

0 


2 7 f. 


-57.55 


1 

0 

0 


Ti\BLS A A COEDFICIS'jiTS H'-I THE: 
BXI2BESI0E (A2.7) 


..LJ 


The location of the oxigni and the 
conic section a,ssociated ?Jith quadrat 
^ and D,. are given in Table A5 . 


inc 1 ina.t i on o f 
c expressions 



D 

JJ^ 

Dp 

Origin of the 
co-iic section 

(a^ s ap) 

-1.753, 

.3620 

-2 .887, 
-.1766 

2.382 , 
-.2715 

Inclination in 
degrees 

-1.322 

3-519 

-3.136 

TABES A5 j 

DETAILS 01 THE 

COLIC SECT] 

:0FS 
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in the transformed coordinates ttlat is the principal 
coordinates of the conic section, S is gi^en as 


D 


0.1617 



3.704 


+ 0.8191 


(A2.8) 


It is seen frctn (A2.8) that D = 0 represents a hyperbola 
Siffiilanly I)^ and are written as 


2^ = - 3.246 I' 


.1 


13.6 + 27.16 


2^ = -.6576 ^^ - 204.2 rf + 7. 811 


(A2.9) 

(A2.10) 


where '^’^7 their respectii^e principal coord ina-tes . 

2-] = 0 and 2p = C therefore represent ellipses, and are 
shown on Figure S. 


2. A short version of the program for the symmetric 
2-terr;i calculation hased on (A2.3) and the coefficients 


..iven in fable 14 > 


be written. 'This is suitable for a 


Quick calculation of quantities such as pressure, skin 
friction and ceritraline velocitv. 



APFEHDIX 3 

2 n£ SECCijD Sn^&uLEii irOII'ix QJ' (;;;). 44) 

L. axLd in equation (3 .44) are given ty (42.5, 6) 

I i 

u-e consider the Dehavioux in thu iieighDournood of the 
second singular point 1. The coordinates of point P are 


■.7309 , agp = -.1849 


(A3.1) 


^'e aeiin? 


>> 

-7 




"IP 


0 up - 


(13.2) 


For the analysis in the neighbourhood of P, the 

are considered, be can then 


linear terms in and 


write 


da. 


G 3 j 




/ O' a> ) ^ -f ( o* D ^ j ^ 

1 * -i' ■ I t di £ 


(43.3) 


no p t p, Y P 


— : i 


1 i ir > 1 


1 /oiSp ^ p £ 




(43.4a) 


(EDo/4a,)^,.p ^ + (4'S„/3a..)^ f 


2 ^ P S2 . 


(43.4b) 


Here d.ash does not denote d ifferent iat ion with respect to r 
but with respect to some other parameter.. Inserting the 
numorical values in the above expressions me get 
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r. = -24.52 

^ I 


-h 2 


5.96 


■>-1 


'1 


- 171 . 4^2 ’ 

+ 33.89 ?2 • 


(13.5 a) 

(13.512) 


‘Ihe above can also le ?iritten as 
.. / 


421 


i -24.52 171.4 1 

{ ; 

i 5.96 33.8si 


I 


■’1 , 
2 ’ 


(13.6) 


Iherefors trace T' and the determinant iTl of the matrix on 
the right side of (13.6) are 


- Q "^7 

— y * J I } 


-1852 . 


/ t ■7. 7 % 


Is .i'KO, P is a saddle point. Phe eigenvalues of the 
above matrix equation are found to be 47.93 and -38.56. 
Ihe principal directions then are, 


‘ ) 1 ~ 

4l - 2.36. 

‘h ^ 

0.082 + 

In t he coord ina t e 

system ( 

written as 


t 

=: ■ 

-38.56 ^). 

/ 1 

72 = 

47.93 


.2 ’ 
'^2 
n ■ 


(13.8a) 

(I3.8t: 


(13. 9a) 


(13.9b) 



APrEHDIX 4 

OK PHE COHPUxSR PROGRAM 


1. The prea icto:c;-co3:rect''.r method Uwsed to integrate 

equation (3.47a) is as follo^i's. The predictor used is the 
Euler's raethod , 



/ a^) H 


(A4.1) 


k'lhere h is the predicted value of a, at the next point 

IP ^ 

of integration, aiid H is the step length in . The 
corrector used is the Euler's modified method, 

s - % = “^2 I (4/4) + (4^4)] ’ 


7 Jhere c^ is the corrected value. The triaximuni (3; ifference 
"between the predicted and the corrected values is given as 


mar^c 

m 



(A4.3) 


In general more than one use of the corrector is 
necessary to get e. within a prescrioed level of accuracy 
ACI". HoweA-er the repeated use of the corrector is permitted 
only upto a prescribed numher of tirr.e i'L4X00E vuthin which if 
e does not fall 'belovi' ACY the step length is halA'ed , E = H/2. 
and the procedure is started again from (A4.1). Honevoi in 


halving the step length if H falls "below a prescribed mini- 
mum ’:.LC'.. , the integration is terminated as the roundoff 


e r r c r w i 11 be s ig n i f i c ant . 
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At the end of the integration over the step, equation 
(3.47h) is integrated as 

Ax = H/2 d/a'^ + l/D^) , (A4.4) 

v'here Ax is the incrernent in x. Since the error of the 

3 

corrector (A4.2) is of order K , the step length is aonhled , 
E = 2H, if e A lCY/10 = ACYIOE , at the end of the integra- 
tion over the step. 

A schematic lajout of the above predictor-corrector 
method nhich incorporates the adjustment of the step length, 
is shown on the next page. At the "beginning of the inte- 
gration a,„ and a axe i-oaowin. Equation (3.46) is used to 

ui m 

determine the independent variable • Ihe diagram explains 
the integration over one step. 

liectaixgular boxes in the diagram showj mathematical 
opexations and the logical operations are shown bv diamonds. 
Print statements are deleted to keep the diagram simple. 
Arrows show the flow of operation, fhe mathematical opera- 
tions are stated within boxes either explicitly or inside 
brackets. Eev' steps of the diagraim are ex planned. 1116 

first box contains b_(l4.l) which means that b are calcula- 

m ffi 

ted using equation (A4.1). To obtain the derivative b' , a 

til ^ 

subroutine YinlPIS pg used, which is shown in another box. 
i. 1 p G X © Gi s. i n operations are clear from the diagram. 
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2. The seq'aerice of operations in a calculation is as 

follov^s; 

Information in the form of v^hether the calculation is 
symmetric or asymmetric, is required to determine the expan- 
sion set. N is the number of terms in (3.33). In case of 
symmetric calculation IT is the number of even eigenfunctions 
included in the calculation. 

The required eigenfunctions and its adjoint functions 
are generated. Coefficients and 12^, given by (3.35) 

are calculated . 

The initial condition for a^,, given by (3.37) is 
c alculated next . 

The integration of the equation (3.36) forms the 
heart of the program, and is carried out in a separate sub- 
routine 'rKASEI'T’. This requires a subroutine YrnllvS to 

T 

solve the set of equations (5.36) for the derivatives a„. 

The predictor-corrector method used for the intef ra,t ion has 
been described above. In abdition follovJing features have 
been incorpoXcireQ . 

i) If the initial condition is such that B <0, EETLW 
is executed with a message. 

ii) If at any stage of the computation D changes sign, 
EEfUiLH is executed with 


a message 
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iii) If maxm. !a^ : is less tjaan or eque.l to a prescritied 
m i ; 

small number ^ , EEi’UEI^ is executed (linear sub- 
range, see Section 3.4.2). 


iv) If the computation proceeds longer then a prescribed 
number of steps, RET UK:! is executed. 

After the successful execution of PHISER, the 
desired results are printed out. 




a. SYMMETRIC SUDDEN EXPANSION 



b. ASYMMETRIC SUDDEN EXPANSION 



c. BASE IN A CHANNEL 


FIG. I TYPES OF SEPARATED FLOW 
PROBLEMS CONSIDERED 



FIG. 2 ORDER SPACE REPRESENTATION OF 
LARGE REYNOLDS NUMBER LIMITS OF NSE 

(Log SCALE FOR L^ , Ly & f AXES) 



FIG.3 DOMAIN OF VALIDITY FOR 
LIMIT EQUATIONS WHEN Ly^l 

(Log SCALE FOR L^ AND Y AXES) 



////,// / // / // // / //,//// / f , f / f / / f ! 



■x/ y / / / / / / V / / / / / / ■/ y!. / / I /. / / e I /■ / / 



FIG. 4 SCHEMATIC OF TWO TYPES OF 
INTERNAL SEPARATED FLOWS 
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FIG. 5 DEPENDANCE OF STREAIV1\MSE 
LENGTH SCALE ON REYNOLDS NUMBER 


B C D 



SEPARATION 

POINT 


FIG. 6 SCHEMATIC OF FORWARD AND 
REVERSED FLOW REGIONS 
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FIG. 8 PHASE PLANE 



















FIG. 15 (Ulc"')^ IN AN INTERMEDIATE 

RANGE OF X 
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FIG-17 REPRESENTATION OF THE ENTRY 
ONDITION FOR SYMMETRIC SUDDEN EXPANSION 


FIG.I8 STREAMLINES FOR SYMMETRIC 
SUDDEN EXPANSION 





FIG-19 CENTRELINE VELOCITY DISTRIBUTION 


FOR SYMMETRIC SUDDEN EXPANSION 



h = 0-5 


PARABOLIC ENTRY PROFILE 
3- TERM 


X= 0 -02 -04 -06 -08 0-1 



FIG. 20 FLOW DEVELOPMENT IN SYMMETRIC 

SUDDEN EXPANSION 





FIG. 23 CENTRELINE VELOCITY DISTRIBUTION 
FOR SYMMETRIC SUDDEN EXPANSION 






FIG. 24 PRESSURE DISTRIBUTION 









FIG. 27 R 





FIG. 28 REPRESENTATION OF THE 
ENTRY PROFILE 



rol- 



FIG. 29 SKIN FRICTION FOR ASYMMETRIC 

SUDDEN EXPANSION 
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FIG. 32 PRESSURE DISTRIBUTION FOR 
ASYMMETRIC SUDDEN EXPANSION 




FIG. 33 REPRESENTATION OF THE ENTRY 
PROFILE FOR THE BASE IN A CHANNEL 
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FIG. 34 REPRESENTATION OF THE 
ENTRY PROFILE 
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